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1.  Project  Goals 

/ 

The  overall  goal  of  this  project  is  to  construct  and  bring  into 

operation  a  free -electron  laser  (FEL)  driven  by  the  racetrack  micro tron  — 

-  r  ) 

(RTM)  electron  accelerator  at  NBS.  NBS  is  providing  the  RTM,  radiation- 
shielded  laboratory  space  for  the  RTM  and  FEL,  and  experimental  areas. 

The  scope  of  the  project  includes:  modifying  the  RTM  injector  for 
increased  peak  current;  developing  electron  beam  transport  from  the  RTM 
to  the  FEL;  developing  a  wiggler  and  optical  cavity;  and  developing 
optical  beam  transport  and  diagnostics.  ^The  plan  for  the  completed 

facility  is  shown  in  Figure  1  of  Appendix  K.  A  jb Wo  '  « 

i  J),tti-»chcv  Jess  ■ 

2.  Activities  and  Accomplishments*  r  ^  rs  •  (-edc  J  ^ — 

In  the  first  year  of  the  project  good  progress  was  made  on  many 
fronts.  The  accomplishments  are  summarized  here  and  are  presented  in 
more  detail  in  the  sections  that  follow.  The  conceptual  design  of  the 
FEL  was  completed,  including  decisions  on  optical  cavity  length  and  pulse 
frequency,  and  detailed  calculations  were  done  to  predict  its 
performance.  For  the  wiggler,  specifications  were  written,  proposals 
from  several  potential  vendors  were  evaluated,  a  vendor  was  selected,  and 
a  contract  was  prepared  for  signing.  The  shielded  room  that  will  house 
the  FEL  was  cleared  of  other  equipment.  Three  schemes  to  increase  the 
RTM  peak  current  were  identified  and  are  being  investigated.  The  5  MeV 
RTM  injector  beam  tests  were  completed,  and  the  electron  transport  line 
from  the  injector  to  the  RTM  was  installed.  The  beam-optical  design  of 


the  electron  transport  line  from  the  RTM  to  the  FEL  was  completed. 
Requirements  for  the  user  facility  were  identified,  and  we  began  clearing 
the  existing  experimental  area  and  designing  the  new  one. 

2.1  System  Design,  FEL  Physics,  and  Performance  Calculations 

On  September  9,  1987  a  conceptual  design  review  of  the  NBS-NRL  FEL 
was  held  at  NBS.  Members  of  the  design  team  presented  all  aspects  of  the 
design  to  the  project  staff  and  to  two  outside  consultants,  Dr.  Richard 
Freeman  of  AT&T  Bell  Laboratories  and  Dr.  Lloyd  Young  of  Los  Alamos 
National  Laboratory.  The  presentations  and  discussions  that  took  place 
are  summarized  in  the  Conceptual  Design  Report,  which  is  attached  to  this 
report  as  Appendix  A.  In  addition,  detailed  design  notes  on  different 
aspects  of  the  FEL  were  written  by  the  project  staff  throughout  the  year. 
These  are  attached  as  Appendices  B-J . 

We  have  used  the  self-consistent,  three-dimensional  computer  code 
SHERA,  which  was  developed  at  NRL,  to  assess  the  effects  of  the 
transverse  emittance  of  the  electron  beam  on  the  performance  of  the  FEL. 

A  complete  description  of  this  work  is  given  in  Appendix  K.  We 
calculated  values  of  the  small-signal  power  gain,  Gp ,  for  normalized 
transverse  emittance,  «n ,  between  5  mmmrad  (or  pm  )  and  20  pm,  with  the 
results  shown  in  Figure  5  of  Appendix  K.  The  design  value  of  en  for  the 
RTM  is  5  pm,  and  the  measured  value  for  the  5  MeV  injector  beam  is  less 
than  1  pm.  The  power  gain  calculated  for  the  design  emittance  will 
readily  support  lasing  over  the  planned  wavelength  range.  It  may  also  be 
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seen  from  Figure  7  of  Appendix  K  that,  due  to  the  good  emittance,  some 
electron-beam  guiding  of  the  optical  beam  can  be  expected  for  the  nominal 
peak  electron  current  of  2  A. 

We  have  developed  a  Gaussian  optical  beam  formalism  for  including 

t 

mirrors  in  SHERA  to  convert  it  to  an  oscillator  code.  This  will  allow  us 
to  study  the  evolution  of  radiation  in  the  optical  cavity  and  evaluate 
performance  of  the  FEL  in  the  presence  of  diffraction  and  mode  mixing. 
This  work  is  described  in  Appendix  L.  Meanwhile,  we  have  used  published 
diffraction  tables  to  estimate  losses  from  the  mirrors  and  intracavity 
vacuum  chamber  apertures.  These  losses  are  negligible  for  UV  and  visible 
wavelengths  and  are  small  for  the  IR  wavelengths  at  which  we  plan  to 
lase.  The  losses  from  diffraction  on  the  vacuum  chambers  of  the  wiggler 
and  the  electron  beam  dump  magnet  are  estimated  to  be  0.6%  and  0.2%  per 
pass,  respectively,  at  10  pm.  These  losses  increase  rapidly  with 
wavelength  and  will  prevent  lasing  at  wavelengths  above  12  pm.  Details 
of  these  calculations  are  given  in  Appendices  G  and  H. 

Imperfections  in  the  wiggler  magnetic  field  can  degrade  FEL 
performance  by  causing  the  electrons  to  "walk"  off  the  radiation  axis  and 
by  altering  the  path  length  of  the  electrons  and  hence  their  phase 
relative  to  the  radiation  field.  In  order  to  evaluate  these  effects,  we 
have  developed  computer  codes  to  model  the  propagation  of  electrons 
through  wiggler  fields  with  errors  and  with  no  radiation  field,  and  have 
checked  the  results  for  agreement  with  analytic  expressions.  We  have 
also  modified  computer  code  SHERA  to  include  wiggler  field  errors.  A 
preliminary  SHERA  calculation  was  performed  for  a  wiggler  with  no 
electron  steering  correction  and  with  0.5%  rms  random  (uncorrelated) 
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field  errors.  The  phase  error  calculated  for  this  case  would  degrade 
performance  significantly.  The  real  wiggler,  however,  will  have  electron 
steering  coils  and  field  errors  that  are  predominantly  correlated,  both 
of  which  will  reduce  phase  errors  substanially.  We  are  currently 
incorporating  these  features  into  the  formalism  and  computer  codes. 

For  a  saturated  FEL  in  the  one -dimensional,  low-gain  model  with  a 
perfect,  untapered,  N-period  wiggler,  the  maximum  efficiency  for 
extracting  power  from  the  electron  beam  is  1/2N.1  Using  this  model,  we 
have  estimated  the  optical  output  power  of  our  FEL  for  commercially 
available  laser  mirrors.  This  work  is  presented  in  Reference  2  and  in 
Appendix  G.  The  calculated  cw  output  power  is  in  the  range  50-200  W  for 
200  nm  £  A  <  7  pm,  and  is  approximately  20  W  for  7  /an  <  A  £  10  pm.  We 
are  performing  more  accurate,  three-dimensional  calculations. 

2.2  Optical  Cavity  and  Wiggler 

The  conceptual  design  of  the  optical  cavity  is  complete.  Details 
are  given  in  Appendix  G.  We  have  chosen  a  two-mirror,  linear  cavity  with 
a  mirror  spacing  of  8.062  m.  This  spacing  gives  a  round  trip  light 
travel  time  of  53.78  ns,  which,  at  the  RTM  injector  frequency  of  74.375 
MHz,  implies  that  there  will  be  two  electron  pulses  and  four  light  pulses 
in  the  optical  cavity  at  one  time.  Since  the  RTM  is  a  cw  accelerator  and 
there  is  gain  in  the  laser  only  when  an  electron  bunch  is  present,  the 
FEL  behaves  as  a  cw,  harmonic-mode- locked  laser.  The  cw  nature  of  the 
accelerator  results  in  very  high  average  output  power  for  the  FEL. 
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For  cavity  mirrors,  we  plan  to  use  multi-layer  dielectric  (MLD) 
coatings  on  transparent  substrates.  An  advantage  of  MLD  coatings  is  the 
ability  to  couple  the  light  out  through  one  mirror  by  using  a  partially 
transmitting/partially  reflecting  MLD  coating.  Compared  to  metal 
mirrors,  the  spectral  bandwidth,  or  tuning  range,  is  narrow.  However,  in 
the  visible  and  UV,  metal  mirrors  cannot  be  constructed  with  the  high 
power  reflection  coefficients  (up  to  R  =  99.9%)  that  are  necessary 
because  ours  is  a  low-gain  laser.  The  cavity  mirrors  must  have  high 
damage  thresholds,  as  we  expect  the  peak  irradiance  to  be  1-2  GW -cm'2, 
and  the  average  power  to  be  200-400  kW-cm*2.  As  discussed  above,  the 
losses  from  diffraction  on  the  cavity  mirrors  and  the  intracavity  vacuum 
chamber  apertures  will  be  negligible  at  ultraviolet  and  visible 
wavelengths  and  small  at  the  infrared  wavelengths  at  which  we  plan  to 
operate . 

The  minimum  wavelength  for  lasing  can  be  reduced  by  lasing  on  the 
third  harmonic  because  Gp  for  the  third  harmonic  exceeds  that  for  the 
fundamental,  as  discussed  in  Appendix  J.  This  is  done  by  maintaining  a 
high-Q  optical  cavity  at  the  desired  wavelength  while  decreasing  the 
electron  beam  energy  by  the  factor  3*%.  We  may  be  able  to  reach  150  nm 
by  this  means. 

The  wiggler  will  have  a  peak  field  of  0.54  T,  130  periods  of  2.8  cm 
for  UV  and  visible  operation,  and  65  periods  for  IR  operation. 
Specifications  for  the  wiggler  were  written  in  June,  1987,  and  a  request 
for  proposals  was  sent  to  potential  wiggler  contractors  in  August. 

Several  proposals  were  received  in  October  and  were  evaluated  by  an 
internal  technical  committee,  with  Dr.  Brian  Kincaid  (AT&T  Bell 
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Laboratories)  acting  as  an  outside  consultant.  Negotiations  with 
prospective  vendors  began  in  November.  A  signed  contract  is  imminent. 

The  room  where  the  wiggler  and  optical  cavity  will  be  located  has 
been  cleared  of  interfering  equipment.  Ve  have  measured  ionizing 
radiation  levels  in  this*  area  during  operation  of  the  NBS  neutron  time- 
of-flight  facility,  when  a  5  ktf  electron  beam  from  the  NBS  110  Mev  linac 
passes  through  the  FEL  area.  The  measured  radiation  levels  are  in  the 
range  of  10-100  lerad  per  year,  which  will  give  the  samarium- cobalt 
wiggler  magnets  an  estimated  useful  lifetime  of  5000  years.  Lifetimes  of 
optical  components  at  this  radiation  level  will  be  1-10  years,  longer 
than  lifetimes  due  to  optical  radiation  damage.  Ve  have  also  initiated 
sound  and  vibration  measurements  in  the  room  in  order  to  determine 
isolation  specifications  for  the  optical  cavity. 

2.3  RTM  and  Electron  Beam  Transport 

The  RTM,  which  is  described  in  Reference  2,  consists  of  a  5  MeV 
injector  feeding  a  microtron.  In  the  microtron,  a  pair  of  180*  end 
magnets  are  used  to  recirculate  the  electron  beam  through  a  12  MeV  rf 
linac  up  to  15  times  for  an  energy  gain  of  up  to  180  MeV.  See  Appendix 
E,  Figure  1.  The  beam  can  be  extracted  from  any  of  14  separate  return 
lines  in  12-MeV  steps.  In  the  past  year,  5  MeV  Injector  beam  tests  were 
completed,  and  the  beam  transport  line  between  the  Injector  and  microtron 
was  installed.  The  12  MeV  linac  was  rf-tested  at  full  power,  and  all 
beam  transport  components  on  the  linac  axis  and  the  end  magnet  vacuum 
chambers  were  installed.  All  magnets  for  the  return  lines  were  designed, 
and  beam-optical  design  of  the  transport  line  between  the  RTM  and  the  FEL 
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was  completed.  We  are  currently  preparing  to  perform  17  MeV  beam  tests 
with  one  pass  through  the  linac. 

The  results  of  the  RTM  Injector  beam  tests  are  given  in  Reference  3. 
The  measured,  normalized  transverse  emittance  at  5  MeV  is  less  than  1  pm, 
considerably  better  than  the  10  pa  required  for  FEL  operation  at 
wavelengths  down  to  90  pa.  The  measured  energy  spread  of  5  keV  at  5  MeV 
is  below  the  criterion  of  0.4%  for  lasing.  Furthermore,  we  expect  the 
fractional  energy  spread  to  decrease  with  increasing  beam  energy. 

The  5  MeV  beam  transport  system  between  the  injector  linac  and  the 
RTM  linac  axis  is  now  complete,  including  alignment  of  all  bending 
magnets  and  quadrupole  focusing  magnets,  and  leak  checks  of  the  vacuum 
envelope.  All  magnets  on  the  accelerator  axis  have  been  installed  and 
aligned. 

The  two  sections  of  the  12  MeV  linac  on  the  RTM  axis  have  each  been 
conditioned  to  a  power  level  between  115  kW  and  120  kW.  This  corresponds 
to  a  voltage  gradient  exceeding  that  required  for  accelerating  the 
recirculating  beam  to  the  nominal  full  energy  of  185  MeV.  A  beam  line 
for  testing  the  beam  after  one  pass  through  the  12  MeV  linac  is  presently 
being  installed  at  the  exit  of  end  magnets  El.  This  beam  line  includes  a 
high-power  beam  stop  and  sufficient  beam  profile  monitors  for  beam 
emittance  and  energy  spread  measurements . 

Designs  have  been  completed  for  the  steerers,  S19  -  S21,  and  for 
quadrupoles  Q10  and  Qtl  on  the  RTM  return  lines.  The  electron  beam 
optics  design  for  the  RTM-FEL  beam  transport  line  is  complete  and  is 
presented  in  Appendix  E.  This  design  includes  proper  transport  of  the 
RTM  beam  over  all  energies  between  29  MeV  and  185  MeV,  and  provides  for 
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an  achromatic  beam  on  the  FEL  axis.  Preliminary  designs  have  been 
established  for  the  movable  deflecting  magnet,  Dll,  and  the  first  two 
dipole  magnets,  D12g  and  D1Zv  The  remaining  magnets  in  the  transport 
line  are  (See  App.  E,  Figs.  1  and  3):  three  2-inch-bore  quadrupole 
doublets  and  a  dipole  bending  magnet  in  the  RTM  room;  two  quadrupole 
doublets  that  act  as  a  telescope  to  form  the  optimal  beam  size  and  waist 
position  on  the  FEL  axis;  an  achromatic  system  of  4  dipoles  to  guide  the 
electron  beam  around  the  upstream  cavity  mirror;  and  a  90°  beam  dump 
magnet  before  the  downstream  mirror. 

2.4  High-Current  Injector 

The  RTM  injector  consists  of  a  100  kV,  5  mA  dc,  thermionic  electron 
gun  followed  by  a  chopper/buncher  system  and  a  5  MeV  cw  linac.  The 
injector  produces  beam  pulses  at  2380  MHz  (the  rf  drive  frequency)  with  a 
maximum  of  0.35  pC  per  pulse.  A  peak  electron  beam  current  of  2-4  A  is 
required  in  the  FEL  to  achieve  lasing.  For  a  3.5-ps  long  beam  pulse, 
this  is  7-14  pC  per  pulse.  We  have  proposed  to  increase  the  charge  per 
pulse  without  increasing  the  average  current  in  the  RTM  (limited  by  the 
rf  power  available)  by  reducing  the  beam  pulse  frequency  to  74.375  MHz, 
the  32nd  subharmonic  of  the  rf  frequency. 

We  are  investigating  three  possible  ways  to  achieve  the  required 
charge  per  pulse:  subharmonic  bunching,  a  laser-driven  photoelectron  gun, 
and  subharmonic  chopping.  In  the  subharmonic  bunching  option,  a  3-ns 
long  pulse  from  a  110  kV  thermionic  gun  goes  through  a  bunching  cavity 
operating  at  74.375  MHz  and  is  compressed  to  70  ps.  The  beam  is  then 
further  compressed  by  the  existing  buncher  to  15  ps  before  entering  the 


Injector  linac,  where  it  is  further  compressed  to  3  ps.  This  technique 
has  the  advantage  of  having  low  source  currents,  but  impresses  a  large 
energy  variation  (±  20%)  on  the  beam  due  to  the  large  bunching  ratios 
needed. 

An  electron  gun  driven  by  a  mode- locked,  cw  laser  requires  a  photo¬ 
cathode  with  a  1%  quantum  efficiency  and  a  one-week  lifetime  to  be 
acceptable  for  user  facility  operation.  We  are  collaborating  with  the 
AT- 7  group  at  Los  Alamos  to  identify  a  suitable  cathode  material. 

Several  materials  have  been  found  with  more-than-adequate  quantum 
efficiency,  but  their  lifetimes  are  short:  a  few  hours  even  in  high 
vacuum  (10-10  Torr)  at  average  currents  that  are  less  than  2%  of  what  we 
require.  Further  tests  at  higher  currents  will  begin  soon  at  Los  Alamos 
with  a  dc  gun  configuration.  Meanwhile,  we  have  been  operating  a  cw, 
mode-locked  Nd : YAG  laser  as  a  prototype  photocathode  driver. 

A  very  promising,  recently-conceived  method  of  increasing  the  peak 
current  involves  replacing  the  present  dc  electron  gun  with  a  100-150  mA, 
110  kV,  thermionic  gun  pulsed  at  74.375  MHz.  The  2 -ns  long  gun  pulses 
would  then  be  chopped  and  bunched  to  the  required  length.  The  required 
electron  gun  is  commercially  available,  the  chopper  would  be  similar  to 
the  existing  one,  and  the  existing  buncher  would  remain  intact.  We  are 
currently  doing  design  calculations  for  this  system.  After  completion  of 
these  calculations  and  the  Los  Alamos  dc  photocathode  tests  early  this 
summer,  we  will  decide  which  method  to  use. 
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2.5  User  Facility 


The  user  facility  will  consist  of  two  experimental  rooms.  The  first, 
XA1 ,  is  an  existing  room  adjacent  to  the  FEL,  about  1600  ft2  in  area  and 
approximately  40  ft  below  ground  level.  It  contained  a  beam  transport 
line  and  a  magnetic  spectrometer  connected  with  the  NBS  110  MeV  linac. 

We  removed  most  of  this  equipment  this  year  for  conversion  to  the  FEL 
facility.  A  building  addition  at  ground  level  will  contain  a  second 
experimental  area,  XA2,  of  about  2400  ft2.  In  the  past  year  we  visited 
several  laser  laboratories  to  define  requirements  for  experimental  areas. 
We  will  take  several  precautions  to  minimize  vibrations  and  sound  in  the 
user  facility.  These  include:  thick  exterior  walls  to  minimize  wind- 
induced  vibrations;  isolation  of  motors,  etc.  from  the  experimental  area; 
and  constant,  laminar-flow  air  conditioning.  We  also  drilled  a  test  hole 
for  location  of  the  optical  transport  system  between  the  two  experimental 
areas.  A  feasibility  study  is  underway  to  determine  if  the  site  provides 
satisfactory  support  for  the  proposed  addition.  A  preliminary  building 
design  was  completed,  and  we  are  preparing  a  bid  package  for  engineering 
and  architectural  design  and  construction. 
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3.  Plans  for  the  Coming  year 


In  the  period  April  1,  1988  through  March  31,  1989,  we  plan  to: 

•  Select  a  method  of  increasing  the  charge  per  pulse  of  the  100  keV 
injector  to  at  least  7  pC. 

•  Design  components  for  a  new  injector  based  on  the  method  selected. 

•  Perform  one-pass  beam  tests  of  the  RTM  linac. 

•  Procure  components  for  the  RTM  return  beam  lines  and  the  RTM-FEL  beam 
transport  line. 

•  Sign  a  contract  with  a  wiggler  vendor  and  oversee  the  magnetic  and 
engineering  design  and  initial  construction  of  the  wiggler. 

•  Calculate  the  PEL  gain  using  computer  program  SHERA  and  a  realistic 
model  of  the  wiggler,  including  steering  and  coherent  field  errors. 

•  Evaluate  FEL  performance  for  a  range  of  electron  pulse  lengths  and 
longitudinal  profiles. 

•  Study  mirror  damage  in  FEL's  and  apply  this  knowledge  to  our  case. 

•  Prepare  an  engineering  design  of  the  optical  cavity. 

•  Prepare  a  conceptual  design  of  optical  beam  transport  to  the  user 
areas . 

•  Complete  the  removal  of  existing  equipment  from  XAl  and  begin 
conversion  to  an  FEL  user  facility. 

•  Complete  the  architectural  and  engineering  design  of  the  building 
addition  and  begin  construction. 
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APPENDIX  A 

NBS-NRL  Free-Electron  Laser  Conceptual  Design  Report 


Injector  Modifications 

Our  goal  is  to  develop  an  injector  with  a  charge  output  of  12  pC  locked 
to  the  32nd  subharmonic  of  the  RF  frequency.  The  desired  beam  parameters 
after  acceleration  in  the  RTM  are  en  (transverse)  <  10  ym,  e  (longitudinal) 

<  30  keV  degrees,  peak  current  4  A,  pulse  length  3  ps.  Frequencies  corres¬ 
ponding  to  the  8th,  16th,  64th  and  128th  RF  subharmonics  are  also  desirable, 
consistent  with  the  average  power  of  the  RTM  klystron.  A  peak  current  after 
acceleration  of  2  A  is  acceptable  for  FEL  performance. 

Our  baseline  design  is  a  Cs3Sb  photocathode  excited  by  a  mode-locked 
frequency -doubled  Nd-Yag  laser.  We  hope  to  make  use  of  the  ongoing  develop¬ 
ment  program  at  Los  Alamos  by  Rich  Sheffield  et  al.  The  Los  Alamos  staff  has 
expressed  a  willingness  to  cooperate  with  us  but  a  formal  collaboration  has 
not  been  established.  The  main  unknown  appears  to  be  photocathode  lifetime 
and  reliability. 

We  recognize  that  the  photocathode  concept  is  in  the  research  and  proto¬ 
typing  stage:  No  working  accelerator  uses  a  photocathode  injector.  A  great 
deal  of  new  information  should  be  available  from  Los  Alamos  in  the  next  few 
months.  If  the  new  information  is  disappointing,  the  possibility  of  a 
thermionic-cathode  gun  with  subharmonic  bunching  appears  to  be  a  viable 
option. 

The  existing  electron  gun  geometry  is  not  useable  at  the  required  current 
level  because  the  anode-cathode  spacing  is  too  large.  Bill  Hermannsfeldt  at 
SLAC  has  calculated  the  performance  of  a  new  gun  design  with  reduced  anode- 
cathode  spacing.  He  predicts  a  normalized  emittance  of  about  2  ym  at  a  cur¬ 
rent  of  0.2  A  at  an  operating  voltage  of  110  keV,  which  would  satisfy  all  of 
our  requirements.  (0.2  A  for  60  ps  at  110  keV  would  compress  to  4  A  for  3  ps 
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after  acceleration.)  The  space -charge  limit  of  Bill's  design  is  about  2  A  at 
110  keV. 


Calculations  must  be  done  with  PARMELA  to  determine  if  the  existing 
capture  section  can  handle  the  higher  current  (by  a  factor  of  40)  and  higher 
gun  voltage.  It  may  be  necessary  to  add  solenoidal  focussing  (B^  field)  over 
the  entire  capture  section.  In  the  worst  case,  the  capture  section  might  have 
to  be  redesigned  and  replaced.  These  calculations  should  be  started  immedi¬ 
ately.  These  calculations  should  include  different  assumed  laser  intensity 
profiles  and  pulse  widths  on  the  cathode.  Uniform  illumination  or  a  gaussian 
profile  may  produce  quite  different  performance. 

High  Peak  Current  Effects 

The  existing  RTM  injector  linac  produces  a  beam  in  which  every  RF  bucket 
contains  the  same  charge,  approximately  0.26  pC.  We  now  want  about  50  times 
as  much  charge  in  the  micropulse,  but  with  only  one  RF  bucket  in  32  filled. 

To  put  the  problem  in  proper  prospective,  several  RF  linacs  have  successfully 
accelerated  micropulses  of  about  10  nC  -  three  orders  of  magnitude  above  our 
goal.  What  is  unique  to  our  situation  is  that  we  are  dealing  with  a  microtron 
instead  of  a  one-pass  linac,  and  our  emittance  and  energy  spread  requirements 
are  very  tight.  The  following  list  enumerates  the  known  potential  effects. 

1.  Wakefields 

2.  Intra-beam  scattering 

3.  RF  heating 

4.  Beam  loading 

5.  Beam  breakup  (BBU) 

6.  Beam-beam  scattering 

7.  Overlapping  micropulsas 


We  now  briefly  consider  these  effects  in  order. 

1.  Wakefields  -  The  longitudinal  wakefield  effect  produces  an  uncompen- 
sable  beam  energy  spread  of  about  1%  in  typical  high  current  S-  and 
L-band  linacs.  IYi  our  case,  with  a  factor  of  103  less  charge  but  an 
aperture  of  a  factor  of  3  smaller  than  other  S  band  accelerators,  we 
expect  a  fractional  energy  spread  of  about  10"4.  Taking  account  of 
the  longitudinal  optics  of  the  RTM,  the  result  of  the  wakefield 
energy  loss  would  be  an  increase  of  the  longitudinal  emittance  from 
an  initial  value  of  20  keV  degrees  to  about  30  keV  degrees  at  185 
MeV.  The  energy  spread  and  phase  spread  would  each  increase  by  about 
20%  in  the  process.  A  more  accurate  calculation  of  both  the  longi¬ 
tudinal  and  transverse  wakefield  effects  is  needed. 

2.  Intra-beam  scattering  -  The  "Touschek  lifetime"  of  our  beam  is  at 
least  three  orders  of  magnitude  longer  than  the  transit  time  of  the 
beam  pulses  from  injector  to  FEL,  so  actual  losses  should  be  negligi¬ 
ble.1  Emittance  growth  due  to  intra-beam  scattering  has  not  been 
calculated  but  we  expect  it  to  be  negligible  since  in  our  estimate  of 
Touschek  lifetime  we  actually  used  the  beam  energy  spread  (10  keV)  as 
the  energy  change  which  would  cause  an  electron  to  be  lost. 

3.  RF  heating  -  This  effect  is  due  to  the  high  frequency  components  of 

the  image  currents  in  beam  tubes.  The  effect  has  been  calculated  for 

the  parameters  of  our  beam  and  wiggler  vacuum  chamber  by  Kelvin  Neil 

of  Lawrence  Berkeley  National  Laboratory.  The  heating  rata  is  a 

trivial  few  milliwatts  per  meter  of  chamber  length. 

1W.  Scharf,  "Particle  Accelerators  and  Their  Uses,"  English  translation  edited 
by  Francis  T.  Cole,  Harwood  Academic  Publishers,  New  York  (1986),  equation 
8.25,  page  540. 


t 

4.  Beam  loading  -  This  depends  primarily  on  the  average  current  because 
the  quality  factor  (Q)  of  the  RF  structure  is  very  large  compared  to 
the  subharmonic  number  of  the  beam  pulses.  The  single-pulse  beam 
loading  is  the  lpngitudlnal  wakefl eld  effect  considered  above.  Beam 
current  fluctuations  on  a  time  scale  which  is  slow  compared  to  Q/f 
but  fast  compared  to  the  RF  system  response  time  (~  0.1  second  due  to 
the  mechanically  driven  power  splitters)  result  in  beam  energy  fluc¬ 
tuations  given  by 

r  ’  '('"'tow  )  r-  ■  (1) 

1  y 

where  Np  is  the  number  of  passes  of  beam  in  the  RTM,  and  I  is  the 
average  beam  current  in  amperes.  In  the  worst  case  (Np  =  15, 

I  =  0.55  mA),  a  1%  Increase  in  current  causes  an  0.18%  decrease  in 
energy.  This  shows  the  importance  of  good  current  stability. 

5.  BBU  -  This  is  also  primarily  an  average  current  effect.  The  8BU 
threshold  for  the  RTM  has  been  estimated  to  be  above  2  mA.  However 
the  subharmonic  structure  has  a  different  Fourier  spectrum  than  a 
beam  with  every  bucket  filled.  One  of  the  harmonics  of  the  beam 
fundamental  frequency  could  be  very  close  to  a  blowup  mode  frequency. 
Then  the  blowup  would  be  driven  by  this  Fourier  component  -  possibly 
much  larger  than  the  driving  term  due  to  noise  on  a  true  CW  beam. 
Further  calculations  should  be  performed. 

6.  Beam-beam  scattering  -  The  beam  micropulses  being  accelerated  from  5 
to  17  MeV  "collide"  with  the  pulses  travelling  in  the  opposite  direc¬ 
tion.  We  calculated  this  effect  for  the  every-bucket-filled  case. 

The  net  effect  was  equivalent  to  a  (defocussing)  focal  length  of  a 
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few  thousand  meters  during  the  5  to  17  MeV  pass  (and  even  weaker  at 
higher  energies).  Since  the  sum  of  the  Impulses  after  32  RF  periods 
depends  only  on  the  average  current,  our  first  guess  is  that  this  is 
negligible.  Further  study  is  recommended  (low  priority). 

7.  Overlapping  micropulses  -  When  operating  at  the  32nd  subharmonic,  a 
snapshot  of  the  beam  pulses  in  the  RTM  linac  section  would  show  that 
of  every  32  RF  buckets  there  are  twenty  empty  buckets,  ten  containing 
one  micropulse  each,  and  two  containing  two  micropulses  (of  different 
energies)  each.  The  only  resulting  effect  that  occurs  to  us  is  an 
increased  space  charge  force,  but  since  this  overlapping  occurs  only 
at  higher  energies,  it  is  negligible,  except  perhaps  as  a  modulation 
of  the  w?kefield-induced  energy  spread. 

In  summary,  while  none  of  the  high  peak  current  effects  appear  serious, 
several  should  be  studied  more  carefully. 


RTM 

The  DoE-funded  RTM  development  will  terminate  in  April  1988,  although  the 
RTM  will  not  be  completed  at  that  time.  It  will  be  completed  with  a  combina¬ 
tion  of  NBS  (STRS)  funding  and  FEL  funding.  The  FEl  funding  will  be  used  for 
the  completion  of  the  electron  beam  line  through  the  wiggler. 

An  outline  of  the  RTM  schedule  as  it  relates  to  the  FEL  program  is  as 
follows: 

December  1987  -  begin  one  pass  beam  tests  (acceleration  to  17  MeV) 

March  1988  -  complete  one  pass  beam  tests,  begin  installation  of  re¬ 

circulation  beam  lines  and  extraction  line 

November  1988  -  begin  RTM  final  beam  tests  (with  existing  thermionic- 
cathode  injector  system) 
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April  1989  -  complete  final  beam  tests,  begin  installation  of  photo¬ 

cathode  injector  system 

January  1990  -  provide  beam  to  FEL  to  begin  FEL  testing 

« 

Electron  Beam  Transport 

The  physics  design  of  the  electron  beam  transport  system  for  the  FEL  has 
been  completed.  The  results  are  reported  in  FEL  Design  Note  #4.*  The  system 
consists  of  the  fallowing  subsystems:  (1)  the  RTM  extraction  line  which 
places  a  non-dispersed  electron  beam  on  the  FEL  axis;  (2)  a  four-dipole  achro¬ 
matic  chicane  to  carry  the  beam  around  the  upstream  optical  cavity  mirror;  (3) 
a  two-quadrupol e-pair  telescope  to  match  the  electron  beam  to  the  optical  mode 
in  the  FEL  cavity;  and  (4)  a  dumping  magnet  which  precedes  the  downstream 
cavity  mirror,  provides  the  capability  of  electron  beam  energy  spectrum 
measurement,  and  conducts  the  beam  to  the  beam  dump  to  dispose  of  it  with 
minimal  background  radiation.  The  transport  system  has  a  normalized  accept¬ 
ance  greater  than  the  63  urn  normalized  acceptance  of  the  RTM.  Conceptual 
designs  exist  for  all  magnets  in  the  system  (8  dipoles  and  5  pairs  of  quadru- 
poles).  Detailed  engineering  design  can  start  as  soon  as'personnel  become 
available. 

FEL  Performance  Calculations 

Three  dimensional  calculations  of  single-pass  power  gain  have  been  per¬ 
formed  for  radiation  wavelengths  in  the  range  x  *  0.15  -  12.0  ym.  The  elec¬ 
tron  beam  energy  range  is  y  -  50  to  y  =  350.  Calculations  were  performed  for 


*FEL  Design  Notes  1-7  are  appended  to  this  report  as  Appendix  A. 


normalized  emlttance  ep  =  5,  10,  and  20  urn  and  peak  currents  of  I  =  0.5,  2, 
and  4  A.  The  results  for  I  =  2  A  are  summarized  in  figure  1.  These  results 
are  in  good  agreement  with  earlier  one-dimensional  estimates.  The  calculated 
gains  are  nearly  independent  of  emlttance  so  that  lasing  can  probably  be 
achieved  down  to  A  *  0.2  urn  even  with  en  as  poor  as  20  pm.  Over  the  entire 
range  of  parameters,  the  gain  is  nearly  linear  in  beam  current.  The  effect  of 
optical  guiding  is  becoming  significant  even  at  I  *  2  A  because  of  the  small 
emittance.  The  effects  of  wiggler  field  errors,  finite  pulse  length,  and  beam 
energy  spread  have  not  yet  been  Included  in  the  calculations,  although  the 
computer  programs  that  have  been  developed  are  capable  of  including  them.  The 
effect  of  the  <  0.1%  energy  spread  of  the  RTM  beam  should  be  negligible.  The 
main  effects  of  the  finite  pulse  length  are  expected  to  be  a  shortening  of  the 
optical  output  pulse  length  and  a  corresponding  decrease  in  energy  per  pulse 
and  time-averaged  power. 

Additional  calculations  are  planned  for  the  duration  of  the  project. 

Some  of  the  topics  to  be  examined  are:  harmonic  yields  (of  spontaneous  co¬ 
herent  light),  efficiency  enhancement  by  wiggler  tapering,  and  lasing  on  the 
third  harmonic. 

Wiggler 

The  wiggler  parameters  are  given  in  Table  I.  Since  pulse  slippage  and 
diffraction  effects  become  serious  for  \  >  2  urn,  we  will  use  the  half-length 
wiggler  above  this  wavelength.  After  a  series  of  visits  and  discussions  with 
FEL  laboratories  and  wiggler  builders,  we  have  concluded  that  a  wiggler  to 
satisfy  our  requirements  can  be  obtained  commercially  via  a  fixed-price  con¬ 
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Figure  1  Power  gain  as  a  function  of  radiation  wavelength  from  the 
3-D  calculations  with  l  =  2  A. 


Table  1  Wiggler  Parameters 


HIGGLER  SPECIFICATIONS 

Wavelength  (  28  mm 

Number  of  Periods:  Full  wiggler  130  +  ends 

Half  wiggler  65  +  ends 

Minimum  Gap  10.0  mm 

Peak  Field  at  Minimum  Gap  >0.54  Tesla 

Permanent  Magnet  Material  (Any)  SmCo  alloy 

Taper:  Adjustable  0  -  0.5  mm/m 

Independently  adjustable,  each  half 

Steering:  Both  ends  and  every  33  periods 

Vacuum  Chamber:  8.4  x  16  mm2  bore,  pairs  of  ports  every 

33  periods 

WIGGLER  TOLERANCES 

RMS  field  error  <  27  gauss  (0.5%  BQ  max) 

Integral  field  errors  (both  planes) 

|  /  Bdz |  <  23  gauss  cm 

33  periods 

Third  harmonic  amplitude  <  0. 1  B 

o 

Transverse  field  variation  <  0.5%  B„  for  I x |  <  5  mm 

o  '  ' 


Optical  Cavity 

The  FEL  optical  cavity  will  have  a  round-trip  transit  time  of  128  RF 
periods.  The  distance  between  cavity  mirrors  is  8.0616  meters.  When  the  full 
length  wiggler  is  used,  the  cavity  con  fi  <;<»'• ’.Li  on  is  symmetric,  so  that  the 
radii  of  curvature  of  the  mirrors  are  equil.  The  cavity  is  asymmetric  when 
the  half-length  wiggler  is  used.  The  baseline  design,  described  in  FEL  Design 
Note  #6,  employs  multi-layer  dielectric  mirrors.  The  downstream  mirror  will 
be  partially  transmitting  (0.5  to  4. 0%)  for  outcoupling.  Dielectric  mirrors 
generally  have  higher  reflectivity  and  damage  thresholds  than  broad-band  metal 
mirrors.  Metal  mirrors  are  too  lossy  to  lase  in  the  UV.  The  disadvantage  is 
that  the  useful  bandwidth  of  dielectric  mirrors  is  less  than  20%,  limiting  the 
on-line  tunability  of  the  FEL.  Several  mirrors  could  be  loaded  into  a  chang¬ 
ing  and  aligning  mechanism  to  permit  a  wider  tuning  range  without  breaking 
vacuum,  but  this  may  be  cumbersome  and  would  certainly  be  expensive. 

Diffraction  losses  in  the  optical  cavity  have  been  estimated,  using 
previously-determined  vacuum  system  apertures.  The  vacuum  chamber  of  the 
wiggler  causes  a  diffraction  loss  of  0.6%  at  10  urn,  which  decreases  rapidly 
with  decreasing  optical  wavelength.  All  other  diffraction  losses  are  negligi¬ 
ble.  The  largest  mirror  size  required  is  <  50  mm  diameter  for  the  downstream 
mirror  at  X  =  10  urn.  Over  most  of  the  operating  range,  25  mm  diameter  mirrors 
will  be  adequate. 

The  output  power,  fluence,  and  irradiance  have  been  calculated  for  the 
above  cavity  design  in  FEL  Design  Note  #6.  The  calculation  uses  reflectivity, 
absorption,  and  scattering  coefficients  of  commercially  available  dielectric 
mirrors.  In  these  calculations,  the  extraction  efficiency  is  taken  to  be 
(2N)"1,  where  N  is  the  number  of  wiggler  periods,  which  is  the  theoretical 
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value  for  an  untapered  wiggler.  Wiggler  magnetic  field  errors  will  tend  to 
decrease  the  efficiency  somewhat,  whereas  tapering  the  wiggler  can  increase 
it  dramatically.  For  the  majority  of  the  output  spectrum,  the  outcoupling 
fraction  is  chosen,  as  a  function  of  wavelength,  to  allow  at  least  a  3%  margin 
between  calculated  small  signal  gain  (one-dimensional  calculation)  and  total 
losses.  This  margin  decreases  to  about  1%  for  0.2  nm  <  X  <  0.3  nm.  Graphs  of 
output  power,  fluence,  and  irradiance  are  give  in  figures  2  and  3. 

The  possibility  of  using  an  on-axis  hole  coupler  in  place  of  a  partially 
transmitting  mirror  will  be  investigated.  The  desired  outcoupling  is  small 
enough  that  a  hole  coupler  would  not  seriously  distort  the  optical  mode  in  the 
cavity.  It  is  technically  feasible  to  laser-drill  a  50  urn  diameter  hole  in  a 
1/4"  thick  quartz  substrate.  This  is  about  5 %  of  the  mirror  spot  size  of  the 
fundamental  (X  =  500  nm).  This  is  the  only  known  outcoupling  method  for  ex¬ 
tracting  harmonic  radiation  below  about  150  nm. 

The  possibility  of  a  cavity  dumped  configuration  for  FEL  operation  was 
discussed  briefly.  The  cavity  dumper  could  consist  of  an  electro-optic 
"polarization  shutter"  or  an  acousto-optic  "beam  deflector."  Rick  Freeman 
pointed  out  that  the  acousto-optic  method  is  superior  in  terms  of  damage 
thresholds  and  wavelength  selection,  although  it  may  mean  that  a  three-mirror 
optical  cavity  is  necessary.  This  area  will  be  investigated  further. 

The  most  serious  concern  in  cavity  development  is  mirror  damage.  Serious 
damage  has  been  observed  at  other  FEls.  Harmonic  radiation  in  the  far  UV 
appears  to  be  a  major  factor.  Rick  Freeman  expressed  an  interest  in  collabo¬ 
rating  on  an  experimental  program  to  measure  UV  damage  to  our  mirrors,  using 
the  FEL  on  the  VUV  ring  at  the  NSLS  facility. 


V.' 
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Figure  2  Predicted  average  output  power  as  a  function  of  wavelength.  The  solid  curves  are  the  upper 
limit,  <Pg>/2N,  for  the  case  of  perfectly  reflecting  cavity  mirrors  and  no  diffraction  losses 

The  circles  (for  the  FLW)  and  diamonds  (for  the  HLW)  are  calculated  values  of  output  power 
using  the  published  reflectivities  of  commercial ly-available  dielectric  mirrors,  realistic 
outcoupling  fractions  (as  discussed  in  the  text),  and  calculated  values  of  D. 
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Figure  3  Predicted  output  fluence  per  micropulse  and  peak  irradiance  as  a  function  of  wavelength, 
corresponding  to  the  average  power  in  figure  2.  The  mode  areas  used  to  obtain  these 
values  are  the  beam  1/e  (amplitude)  sizes,  corrected  for  the  defocussing  of  the  output 
coupler.  • 


Design  of  the  optical  beam  transport  system  will  be  deferred  until  we 
have  better  information  on  user  requirements.  Becauseof  the  long  distances 
involved  (up  to  ~  100  m),  the  full  power  beam  will  be  transported  in  vacuum. 
This  will  also  help  with  the  very  severe  safety  concerns  due  to  the  high 
power.  A  safety  protocol  will  have  to  be  developed. 

Optical  Diagnostics 

A  wide  variety  of  optical  instrumentation  is  needed  because  of  the  large 
wavelength  interval,  high  power,  and  short  time  structure.  A  very  wide  dynam¬ 
ic  range  in  intensities  must  be  covered  since  we  must  be  able  to  diagnose  the 
spontaneous  emission  from  the  FEL  in  order  to  tune  up  for  lasing.  For  pulse 
lengths  of  one  picosecond  or  less,  autocorrelation  methods  appear  to  be  the 
only  available  means  for  temporal  measurements. 

Electron  Beam  Diagnostics 

Methods  for  measuring  electron  beam  size  profile  and  emittance  are  avail¬ 
able.2  The  ability  to  measure  the  beam  energy  distribution  will  be  incorpor¬ 
ated  in  the  beam  dump  magnet.  A  method  for  fast  temporal  measurements  is  not 
readily  available.  Streak  cameras  have  time  resolutions  of  a  few  picoseconds, 
which  is  not  adequate  for  optimizing  the  peak  current.  The  possibility  of 
using  transition  radiation  interferograms  for  temporal  measurements  will  be 
investigated. 


2M.A.  Wilson,  et  £^.,  "Performance  of  the  5  MeV  Injector  for  the  NBS-Los 
Alamos  Racetrack  l-licrotron,"  1987  IEEE  Particle  Accelerator  Conference,  IEEE 
Catalog  37CH2387-9,  (1987)  pp  322-324. 
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Control s 

The  control  system  for  the  FEL  will  be  integrated  with  the  existing  RTM 
controls.  Unfortunately,  the  RTM  control  system  is  rapidly  approaching  its 
maximum  capacity.  It  is  not  feasible  to  add  a  dedicated  secondary  station  for 
the  FEL.  Both  hardware  and  software  improvements  are  needed.  In  the  hardware 
area,  the  major  change  needed  is  to  convert  the  8-bit  microprocessors  to  32- 
bit.  In  addition  to  increasing  the  processing  speed,  this  will  allow  us  to 
expand  the  data  base  beyond  the  present  limit  of  1200  devices.  This  change 
would  require  a  complete  rewrite  of  the  secondary  station  software  and  the 
redesign  and  fabrication  of  all  special  interface  boards,  an  estimated  3-5  man 
year  effort.  A  gradual,  step-by-step  upgrading  of  the  system  is  planned  over 
several  years.  This  will  begin  with  the  purchase  of  a  32 -bit  dual  processor 
system  for  the  yet  to  be  constructed  primary  station  wire  scanner  display 
subsystem.  This  32-bit  dual  processor  system  will  not  only  form  the  wire 
scanner  display  system,  but  has  sufficient  processing  power  to  permit  the 
replacement  of  the  present  control  system  DEC  POP  11/44  CPU  and  8-bit  CP/M 
system.  This  can  be  done  with  a  minimum  amount  of  effort  (~  1  man  year). 

Later  addition  of  a  third  processor  to  this  system  would  complete  the  replace¬ 
ment  of  the  existing  8-bit  processors  in  the  primary  station.  This  would 
require  a  much  larger  effort  (~  2-3  man  years)  to  recode  the  existing  primary 
station  hardware  and  replace  the  primary-to-secondary  station  link.  It  is 
essential  that  we  have  a  functional  system  for  beam  tests  in  November  1988  and 
in  January  1990. 
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Building 

The  existing  above  ground  building  that  houses  the  neutron  time  of  flight 
facility  (NTOF)  will  be  modified  with  a  2000  ft2  addition  for  FEL  u.sers. 

Except  for  a  relatively  small  office  area,  the  addition  will  be  configured  as 
a  single  large  research  hall.  Partitions,  optical  curtains,  etc,  will  be 
installed  to  suit  user  needs.  Staff  visits  to  other  laboratories  will  help 
us  plan  for  user  needs.  Detailed  engineering  design  of  the  addition  is  being 
carried  out  by  the  NBS  Plant  Division. 
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FEL  Design  Note  #1 


Issues  Involved  in  the  Choice  of  Injector  Pulse  Frequency 
and  Optical  Cavity  Length 

S.  Penner 

.  April  30,  1987 

The  round  trip  transit  time  must  be  an  integral  multiple  of  the  RTM  RF 
period.  Thus  the  cavity  length  L  is  given  by 

L  -  N  c/2f0  (1) 

where  f0  *  2.38  x  109  Hz,  C  =  2.9979  x  108  m/s,  and  N  is  any  integer.  For  N  = 
120  (L  *  7.558  m)  we  have  determined  that  the  optical  cavity  will  fit  in  the 
available  space.  It  will  be  difficult  to  work  with  a  shorter  cavity.  A  some¬ 
what  larger  cavity  is  feasible,  especially  if  the  cavity  mirrors  can  be 
asymetical  with  respect  to  the  wiggler  center,  or  if  we  can  reduce  the  shield¬ 
ing  wall  (between  RTM  room  and  magnet  room)  to  six  feet.  A  longer  cavity 
reduces  the  heat  load  on  the  mirrors,  but  makes  the  mirror  alignment  more 
difficult.  N  =  128  (L  *  8.062  m)  appears  to  be  feasible. 

The  injector  pulse  frequency  fj  must  be  an  integral  sub-multiple  of  the 
RF  frequency  and  an  integral  multiple  of  N;  i .e.  fj  =  m  f0/N,  where  m  =  1,  2, 
....  We  want  to  be  able  to  operate  at  m  »  1,  which  corresponds  to  the  minimum 
average  beam  power,  and  at  the  largest  m  that  is  compatible  with  the  available 
electron  beam  power,  P  =>  100  kW.  For  our  minimum  peak  current  and  micropulse 
length  requirements,  I  *  2A  and  t  =  3  pS,  and  design  maximum  energy,  E  * 

185  MeV, 

P  -  ($)  <  100  kW.  (2) 

Thus  ^  >  26.4. 

m  — 
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Our  original  choice  was  m(raax)  ■  5,  with  N  -  120,  but  it  may  prove  difficult 
to  operate  the  drive  laser  for  the  photocathode  at  m  *  5.  It  may  in  fact  be 
much  easier  and  operationally  convenient  to  provide  m  -  1 ,  2,  and  4  by 
acousto-optical  techniques  applied  to  the  mode-locked  drive  laser. 

The  frequency  ratio  N/m  must  have  excellent  phase  stability,  approximate¬ 
ly  1°  of  f0.  This  may  be  easier  to  achieve  at  some  particular  values  of  N/m. 
Possible  choices  Include,  N/m  *  32  *  25  and  N/m  *  30  *  5  x  3  x  2.  Phase  sta¬ 
bility  depends  not  only  on  the  choice  of  N/m  but  on  the  choice  of  the  master 
oscillator  frequency,  which  could  be  f0,  using  frequency  division  to  obtain 
f.;  or  fj,  using  frequency  multiplication  to  obtain  f0.  Finally,  the  mode- 
locked  laser  frequency  should  be  in  the  range  where  commercial  Instruments  are 
available,  roughly  80  -  140  MHz. 

Some  numerical  possibilities  are: 


N 

m(max) 

N/m(max) 

Other  Possible  m  values 

120 

5 

24 

99.167 

1  (?) 

120 

4 

30 

79.333 

1.  2 

128 

8 

16 

148.750 

1,  2,  4 

128 

4 

32 

74.375 

1,  2 
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Our  original  choice  was  m(raax)  *  5,  with  N  =  120,  but  it  may  prove  difficult 
to  operate  the  drive  laser  for  the  photocathode  at  m  »  5,  It  may  in  fact  be 
much  easier  and  operationally  convenient  to  provide  m  *  1 ,  2,  and  4  by 
acousto-optlcal  techniques  applied  to  the  mode-locked  drive  laser. 

The  frequency  ratio  N/m  must  have  excellent  phase  stability,  approximate¬ 
ly  1°  of  f0.  This  may  be  easier  to  achieve  at  some  particular  values  of  N/m. 
Possible  choices  include,  N/m  *  32  *  25  and  N/m  *  30  *  5  x  3  x  2.  Phase  sta¬ 
bility  depends  not  only  on  the  choice  of  N/m  but  on  the  choice  of  the  master 
oscillator  frequency,  which  could  be  f0 ,  using  frequency  division  to  obtain 
f j ;  or  f.,  using  frequency  multiplication  to  obtain  f0.  Finally,  the  mode- 
locked  laser  frequency  should  be  in  the  range  where  commercial  Instruments  are 
available,  roughly  80  -  140  MHz. 

Some  numerical  possibilities  are: 
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RTM  Shielding  Requirements 
FEL  Design  Note  #2 

S.  Penner 
May  12,  1987 

» 

Because  space  in  the  magnet  room  for  the  FEL  Is  very  limited,  the  ques¬ 
tion  has  been  raised  if  we  can  reduce  the  thickness  of  the  shielding  wall 
between  the  RTM  room  and  the  magnet  room  from  12  ft  to  6  ft.  The  original 
design  of  the  llnac  complex  called  for  shielding  sufficient  to  reduce  the 
exposure  in  occupied  areas  to  less  than  2.5  mrem  per  hour  under  "worst  condi¬ 
tions."  Worst  conditions  were  taken  to  be  100  kW  of  110  MeV  electrons  hitting 
an  unshielded,  thick,  high-Z  target.  This  turned  out  to  require  a  shield  of 
12  ft  of  concrete,  assuming  a  target-to-exposure-site  distance  of  24  ft.  The 
calculations  are  contained  in  my  old  data  book  "Design  of  Magnet  System  for 
Linac,"  No.  3,  March  l  -  September  6,  1960,  pages  74-83.  These  old  calcula¬ 
tions  are  experimentally  verified  in  the  sense  that  measured  neutron  fluxes  in 
the  linac  complex  due  to  shield  penetrations  are  below  the  predicted  level, 
often  by  a  large  factor.  One  large  factor  of  conservatism  is  the  assumption 
that  the  fast  neutrons  which  survive  thick  shields  are  isotropic.  In  fact, 
they  are  strongly  peaked  in  the  electron  bean  forward  direction.  For  both 
background  calculations  discussed  here,  the  direction  of  interest  is  >  90° 
from  the  electron  beam  direction,  which  should  provide  a  substantial  addition¬ 
al  decrease  in  dose. 

Now,  we  are  concerned  with  exposure  in  the  RTM  room  when  the  linac  is 
operating.  The  linac  now  operates  <  1200  hours/year  at  <  5  kW  beam  power  at 
energies  <  110  MeV.  The  beam  line  for  neutron  time-of-fl ight  work  Is  36  ft 
from  the  nearest  occupyable  point  in  the  RTM  room.  My  old  calculations  con¬ 
clude  that  the  neutron  flux  through  a  6  ft  shield  due  to  a  110  MeV  electron 
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beam  Is  1.8  x  106  neutrons/kw  sec.  Thus  the  flux  due  to  5  kW  of  beam  at  36  ft 
from  the  source  Is  0.6  neutrons/cm2  sec  *  0.15  mrem/hr.  Since  the  required 
maximum  dose  rate  for  40  hour/week  occupancy  remains  2.5  mrem/hr,  a  6  ft 
shield  is  safe  by  a  factor  of  17.  In  terms  of  Integrated  exposure,  sln.e  the 
llnac  operates  less  than  1200  hours  per  year,  and  we  cannot  expect  any  Indi¬ 
vidual  to  be  in  the  area  more  than  20%  of  his  working  time,  and  the  beam  loss 
in  the  magnet  room  is  <  20%  of  the  total  beam  power,  the  expected  maximum 
exposure  per  year  is  7  mrem.  This  is  a  factor  of  ~  400  below  legal  limits. 

We  conclude  that  reducing  the  shielding  wall  thickness  to  6  ft  Is  completely 
safe. 

We  can  also  calculate  the  exposure  In  MR#2,  which  will  be  the  new  "UV 
Laboratory"  due  to  the  100  kW  RTM  beam  being  dumped  In  the  magnet  room.  The 
minimum  shield  thickness  is  equivalent  to  12  ft  of  concrete  (It  Is  actually 
4  ft  of  concrete  plus  12  ft  of  fill  dirt),  and  the  source-to-exposure  distance 
is  >  24  ft.  Assuming  200  MeV  operation,  the  neutron  flux  Is 
4.3  x  106  n/kw-sec  implying  64  neutrons/cm2  sec,  equivalent  to  160  mrem/hr. 

The  existing  beam  dump  provides  a  reduction  of  a  factor  of  <  100  due  to  the 
fact  that  the  "target"  is  water  rather  than  a  hlgh-Z  material  and  there  Is  a 
4  inch  Pb  shield  built  Into  the  dump.  This  reduces  the  maximum  exposure  In 
one  corner  of  the  UV  laboratory  to  1.6  mrem/hour.  It  will  be  necessary  to 
survey  the  UV  lab  with  the  RTM  running.  Some  additional  shielding  around  the 
beam  dump  may  be  needed.  (Note  that  the  half-value  thickness  of  Pb  for  high 
energy  neutrons  Is  about  1.7  Inches). 
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FEL  Design  Note  #3 
Optical  Cavity  Length  Choice 


S.  Penner,  C.  Johnson,  and  P.  Debenham 
June  19,  1987 

» 

Further  study  of  the  Issues  discussed  In  FEL  Design  Note  #1  have  led  to 

the  choice  of  an  optical  cavity  whose  round  trip  length  in  units  of  RF 

★ 

wavelengths  is  M  *  128,  Instead  of  M  *  120  as  stated  in  the  proposal.  The 
mode-locked  laser  that  drives  the  photocathode  will  operate  at  n  =  32,  i.e., 
the  32nd  subharmonic  of  the  RF  frequency.  Instead  of  the  24th  subharmonic  as 
originally  proposed.  This  choice  corresponds  to  the  bottom  line  of  the  table 
in  FEL  Design  Note  #1. 

Reasons  for  these  choices  are: 

1.  The  longer  optical  cavity,  L  *  8.06158  m,  allows  a  36  inch  working 
space  between  the  back  of  the  dumping  magnet  and  the  downstream 
cavity  mirror.  The  increased  distance  is  needed  for  optical 
diagnostic  and  alignment  equipment,  and  possibly  for  outcoupllng 
and/or  a  device  to  permit  Q-switched  operation  of  the  FEL. 

2.  The  longer  cavity  allows  more  space  around  the  upstream  cavity  mirror 
for  the  electron  beam  chicane.  In  particular,  the  transverse 
distance  between  the  optical  axis  and  the  electron  beam  center  line. 
In  the  plane  of  mirror,  can  be  increased  to  2.97  inches,  compared  to 
1.87  Inches  In  the  old  design.  There  Is  16.2  inches  between  the 
upstream  mirror  and  the  position  where  the  transverse  distance 
between  the  electron  beam  and  the  optic  axis  Is  1  Inch,  so  that 
optical  components  could  be  installed  in  this  region  as  well. 

*We  have  changed  notation  from  Design  Note  #1,  switching  from  N  to  M  for  the 
optical  cavity  round  trip  length,  and  defining  n  as  the  ratio  of  RF  frequency 
to  mode-locked  laser  frequency. 
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3.  The  lower  photocathode  frequency,  74,375  ftiz  (Instead  of  99.167  MHz 
for  n  *  24),  allows  a  peak  current  of  2.42  A  in  a  micropulse  length 
of  3  pS  at  the  electron  beam  design  power  limit  of  100  kW  at  185  MeV 
maximum  energy.  The  older  design  would  have  allowed  only  1.62  A  with 

t 

all  other  parameters  unchanged.  The  higher  peak  current  is  very 

Important  at  high  energies  In  order  to  obtain  the  gain  needed  for 

lasing  at  short  wavelength.  We  have  determined  that  a  mode-locked 

laser  frequency  of  74.375  MHz  Is  not  a  problem. 

★ 

4.  The  m  3  4  design  now  used  allows  for  operation  of  the  photocathode 
with  m  3  2  or  1  by  dividing  the  mode-locked  laser  frequency  by  powers 
of  2.  This  can  be  done  by  electro-optical  techniques  which  maintain 
good  timing  stability.  The  old  design  was  for  m  =  5.  The  only  other 
synchronous  value  would  have  been  m  3  1,  but  dividing  the  optical 
frequency  by  5  is  difficult.  It  Is  very  valuable  to  have  the 
flexibility  of  micropulse  frequency  so  that  the  FEL  can  be  turned  on 
and  tuned  up  with  reduced  electron  beam  power.  Also,  if  we  succeed 
in  obtaining  peak  beam  current  higher  than  2.4  A,  we  can  reduce  the 
micropulse  frequency  to  stay  within  the  beam  power  limit.  For 
example,  at  m  3  2  (37.1875  Wz  micropulse  frequency),  I  3  4.8  A  in 

3  pS  pulses  at  185  MeV  corresponds  to  100  kU. 

5.  Space  in  the  Magnet  Room  for  the  longer  cavity  can  be  provided  by 
reducing  the  thickness  of  the  shielding  wall  between  the  RTM  room  and 
the  Magnet  Room  to  6  ft.  FEL  Design  Note  #2  considers  the  radiation 
safety  consequences  of  this,  and  concludes  that  It  Is  safe. 

*m  Is  the  number  of  electron  beam  micropulses  per  optical  cavity  round  trip 
length. 
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Figure  1  shows  the  accessible  electron  beam  peak  current  that  can  be 
achieved  as  a  function  of  electron  beam  energy  for  various  values  of  m.  The 
upper  energy  limit  Is  set  by  the  magnetic  field  capability  of  RTM  components. 
The  Injector  limit  of  4  A  Is  our  design  goal  for  the  new  photocathode 
Injector.  The  hyperbolic 'curves  labeled  m  *  4,  8,  16  are  Imposed  by  the 
100  kU  limit  on  electron  beam  average  power.  (Values  of  m  >  4  can  be  obtained 
by  splitting,  delaying,  and  recombining  the  optical  pulses  from  the  mode 
locked  laser.)  For  m  <  4,  beam  power  >  100  kW  cannot  be  achieved  within  the 
design  range  of  the  energy  and  peak  current.  For  m  >  16  another  potential 
limit  at  low  energies  Is  due  to  average  current.  The  2  mA  average  current 
value  Is  a  guess  of  the  beam  breakup  (BBU)  limit.  Another  potential  average 
current  limit  might  be  the  average  power  available  from  the  mode-locked 
laser. 

The  power  limit  curves  In  figure  1  are  obtained  from  the  relation 


I 


M  P 

3  m  f0  t  E 


(1) 


where  I  is  the  peak  current,  P  the  power  and  E  the  energy  of  the  electron 
beam,  fQ  is  the  RF  frequency  (2380  MHz),  and  t  is  the  micropulse  length 
(~  3  pS).  The  dependence  of  peak  current  on  average  current,  1,  Is 


I 


M  1 

t  m  f0 


(2) 


Figure  2  shows  the  new  layout  of  the  upstream  optical  cavity  mirror  and 
the  dipole  magnets  near  the  mirror.  The  dipole  magnets  are  made  smaller  than 
our  previous  estimated  size  by  using  water-cooled  colls,  and  a  picture  frame 
design.  There  will  be  a  total  of  5  magnets  of  this  design:  013  which 
completes  the  RTM  extraction  beam  line;  D14-D15,  located  In  the  RTM  room;  and 
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D16-D17,  shown  In  figure  2.  The  four  magnets  014,  15,  16,  17  form  an 
achromatic  chicane  whose  purpose  Is  to  deflect  the  electron  beam  around  the 
cavity  mirror.  The  smaller  transverse  dimension  of  016  results  In  a  clear 
llne-of -sight  along  the  FEL  optic  axis. 

Figure  3  shows  the  new  layout  of  the  FEL  In  the  magnet  room.  This  layout 
Incorporates  the  8  meter  optical  cavity,  the  reduced  shielding  wall  thickness 
and  the  new  design  of  the  achromatic  chicane. 
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FEL  Design  Note  #4 
Electron  Beam  Transport  Design 

S.  Penner  and  E.  Lind strom 
•  July  7,  1987 


I.  Introduction 

This  note  presents  the  physical  layout  and  electron  beam  optics  design 
for  the  RTM/FEL  beam  transport  system.  The  function  of  this  system  Is  to 
transport  the  beam  from  the  RTM  through  the  FEL  wlggler  and  Into  the  electron 
beam  dump.  The  system  consists  of  two  major  sections,  the  RTM  extraction  line 
and  the  FEL  transport  section. 

The  RTM  extraction  line  accepts  the  beam  which  Is  deflected  from  the  RTM 
by  the  moveable  extraction  magnet.  Dll,  and  transports  It  onto  the  FEL  axis. 
The  three  major  criteria  for  this  beam  line  are:  (1)  It  must  transport  the 
beam  through  and  around  a  maze  of  other  RTM  components  (especially  the  RF 
drive  system),  (2)  it  must  function  properly  at  all  beam  energies,  29-185  MeV, 
and  (3)  it  must  provide  an  achromatic  beam  on  the  FEL  axis.  A  design  which 
meets  these  criteria  Is  presented  In  section  III  of  this  note. 

The  FEL  transport  section  operates  on  the  beam  which  exits  the  RTM 
extraction  line.  The  beam  must  be  guided  around  the  FEL  upstream  optical 
cavity  mirror  achromatically,  transported  through  the  wlggler  with  a  beam 
envelope  which  is  concentric  with  the  optical  mode  shape,  and  removed  to  a 
beam  dump,  avoiding  the  downstream  optical  cavity  mirror.  In  addition  to  the 
requirement  of  operating  over  the  full  range  of  beam  energies,  this  section 
must  function  properly  over  the  full  range  of  wlggler  field  strength,  and  must 
provide  the  diagnostic  capability  of  measuring  the  electron  beam  spectrum  with 
and  without  lasing.  The  Resign  of  the  FEL  transport  section  Is  presented  in 
section  IV  of  this  note. 
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II.  Electron  Beam  Model 


In  order  to  calculate  beam  sizes,  aperture  requirements,  and  focusing 
strengths  we  need  a  characterization  of  the  six-dimensional  phase  space  of  the 
electron  beam  In  the  RTM.  ,  At  present,  the  only  available  experimental  data  Is 
the  beam  emlttance  measured  at  5  MeV.1  The  normalized  transverse  emlttance  is 
0.56  t  0.1  mm  mrad  in  the  horizontal  plane,  and  0.70  ±  0.17  mm  mrad  in  the 
vertical  plane.  The  longitudinal  emlttance  has  an  upper  limit  of  5  keV 
degrees.  These  are  essentially  envelope  emlttances.  containing  >  95%  of  the 
beam  current.  However,  we  expect  larger  emlttances  from  the  completed  RTM 
because  of  (at  least)  two  Important  effects.  First,  the  high  current  Injector 
needed  for  the  FEL  will  have  larger  emlttance  than  the  present  Injector. 
Second,  the  Injection  transport  system  from  the  5  MeV  llnac  to  the  RTM  has 
significant  chromatic  aberrations.  For  purposes  of  estimating  beam  sizes  and 
focal  strengths,  we  will  continue  to  use  our  design  goal  emlttances  of 
eN  *  5  pm  transverse  (normalized)  and  20  keV  degrees  longitudinal. 

For  the  purpose  of  determining  the  aperture  requirements  of  the  system, 

we  must  be  much  more  conservative.  The  limiting  aperture  of  the  RTM  is  the 

8  m  llnac.  Its  aperture  Is  5  mm  radius.  With  a  5  MeV  injected  beam  and 

12  MeV  energy  gain  in  the  first  pass,  the  linac  acceptance  Is  e  *  63.4  um, 

m«x 

normalized.2  Our  design  goal  for  the  normalized  acceptance,  AN,  of  the  trans¬ 
port  system  is 

An  >  e  a  63.4  um,  (1) 

n  max 

so  that  the  system  will  never  Intercept  any  beam,  assuming  that  there  are  no 
significant  sources  of  emlttance  growth  after  the  first  pass  through  the 
llnac.  All  beam  transport  calculations  are  performed  with  a  beam  emlttance  of 
eN  *  5  um,  and  the  acceptance  calculated  from 

An  -  eN  •  smallest  (RA/»*C)2,  (2) 
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where  r£  are  the  calculated  bean  envelope  radii,  and  Rft  the  design  aperture 
radii  at  locations  throughout  the  system  In  both  transverse  planes.  From 
equations  (1)  and  (2)  the  minimum  aperture  radii  required  are  given  by 


*min  > 


•  rc  ■  3.56  rc 


Additional  margin  above  this  value  Is  desirable  to  allow  for  misalignments  and 
steering  errors. 

The  orientation  of  the  phase  ellipsoid  at  the  entrance  to  the  extraction 

system  Is  needed  to  calculate  beam  sizes.  The  basic  optics  design  of  the  RTM 

produces  an  upright  orientation  of  the  phase  ellipses  at  the  center  of  the 

accelerating  section  on  every  pass  In  all  three  phase  planes.  We  use  this  as 

the  starting  point  of  our  calculations,  with  the  beam  sizes  In  all  four 

transverse  coordinates  (x,x',y,ya)  taken  from  calculations  performed  with  the 

(first-order)  optics  code  RTM.OPTICS  for  transverse  tunes  of  u  *  u  *  45°  per 

x  y 

pass.3  The  longitudinal  phase  ellipse  coordinates  are  taken  from  STRACE 
calculations  which  give  $  *  1.11  degrees,  6E  ■  18  keV  for  an  area  of  20  keV 
degrees.'*  From  this  starting  point,  the  following  transformations  are  applied 
to  the  beam: 

(1)  The  transverse  beam  coordinates  are  multiplied  by  (1  -  1  2t 

where  p  Is  the  beam  momentum  at  extraction,  to  account  for  the 
transverse  optics  of  the  llnac. 

(2)  The  beam  matrix  Is  given  a  transformation  which  differs  from  the 
unit  matrix  only  In  the  RS6  element,  given  by 

3  6.,MeV  sin  r  1()2  gercent  (4) 

k56  p  x  480j  iU  aegree  *  ’ 

which  accounts  for  the  energy  focusing  of  the  llnac.  (9°  Is  the 
resonant  phase  for  u  ■  90°  per  pass). 
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(3)  The  beam  drifts  5.9  meters. 

These  transformations  project  the  beam  to  the  entrance  of  end  magnet  E2.  From 
this  point  onward,  the  system  Is  represented  by  standard  TRANSPORT  elements. 

The  origin  of  floor  layout  coordinates  In  the  TRANSPORT  calculations  Is 
the  point  where  the  llnac  axis  Intersects  the  effective  edge  of  El.  The 
Z-ax1$  points  north  along  the  llnac  axis,  positive  X  Is  to  the  west,  and  posi- 
tlve  Y  Is  up.  The  X  coordinate  of  the  starting  point  Is  adjusted  for  the 
closed  orbit  offset  correction  of  the  Injection  chicane.  Except  for  this 
correction,  the  Injection  chicane  and  Injection  quads  (Q6-Q7)  are  Ignored. 

End  magnets  El  and  E2,  the  extraction  magnet  Dll,  and  D8  are  explicitly 
Included  In  all  calculations.  El  and  E2  are  constrained  to  have  no  vertical 
focusing  for  normal  entry  and  exit  by  setting  their  gap  height  to  zero. 
TRANSPORT  calculates  the  vertical  focusing  for  non-normal  entry  and  exit  of 
El.  All  other  dipoles  have  vertical  focusing  effects  appropriate  to  their  gap 
dimensions  and  fringe-field  type.  (We  distinguish  between  pole-tip  type  and 
picture  frames  by  the  parameter  In  TRANSPORT). 

Table  I  gives  floor  layout  coordinates,  axis  direction  angles,  magnetic 
fields  and  bend  angles,  for  extraction  at  29  MeV  and  185  MeV. 

In  the  185  MeV  calculations,  the  return  leg  quadruples  are  turned  off 
because  they  have  an  unrealistic  large  effect  on  the  dispersion  In  the  beam 
model  used.  In  actual  operation,  the  dispersion  on  the  llnac  axis  is  not  zero 
on  every  pass  because  of  the  return  leg  quadrupoles.  The  return  leg  quadru¬ 
ples  on  the  extraction  orbit  would  be  adjusted  to  match  the  recirculation 
tune.  This  would  result  in  a  dispersion  on  the  extraction  line  close  to  that 
used  In  the  model  with  return  leg  quads  turned  off,  and  a  converging  beam  into 
the  extraction  line.  In  our  calculations,  we  have  varied  Q10-Q11  strengths 
over  a  realistic  range  of  values,  and  determined  that  achromatlclty  and 
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adequately  small  beam  size  can  be  obtained  by  adjustment  of  the  extraction 
line  quadrupoles.  In  the  29  MeV  case,  there  Is  no  constraint  on  the  setting 
of  Q10-Q11,  and  Me  adjust  them  In  the  calculations  to  obtain  the  desired  beam 
sizes.  t 

III.  RTM  Extraction  Line 

This  system  transports  the  beam  from  the  RTM  onto  the  FEL  axis.  Since  we 
want  to  accomplish  this  without  resorting  to  physically  moving  any  magnet 
other  than  Oil,  the  layout  shown  In  figure  1  was  chosen.  The  extraction  pro¬ 
cess  begins  with  a  deflection  of  the  beam  on  the  appropriate  return  leg, 
through  approximately  2°.  After  passing  through  El  and  D8,  the  beam  has 
attained  enough  horizontal  displacement  to  clear  07  on  the  llnac  axis.  At 
this  point,  the  beam  encounters  D12H,  which  completes  the  deflection  Into  the 
extraction  channel.  Before  012H,  the  bean  axis  depends  on  the  extraction 
energy.  After  D12H,  the  beam  axis  Is  Independent  of  which  return  leg  Dll  is 
located  on.  The  beam  Is  then  deflected  vertically  by  D12V,  after  which  it 
traverses  most  of  the  RTM  room  at  angles  8  -  -6.91°,  0  ■  1.88°  relative  to 

a  Jr 

the  Z  axis.  The  final  extraction  line  deflecting  magnet,  D13,  is  located  at 
the  intersection  point  of  this  angled  beam  line  with  the  FEL  axis.  The  FEL 
axis  Is  parallel  to  the  Z  axis  at  x  *  -1.057  m,  y  ■  0.254  m.  In  order  to 
provide  the  required  bend  In  both  planes,  D13  is  rotated  about  the  FEL  axis  by 
-15.26°  from  the  horizontal  plane.  As  can  be  seen  In  figure  1,  this  beam  line 
passes  very  close  to  some  RF  plumbing  components  of  the  RTM.  Some  adjustment 
at  assembly  time  may  be  needed.  The  location  of  the  FEL  axis  and  the  direc¬ 
tion  of  the  angled  beam  can  be  changed  slightly.  If  necessary.  After  these 
calculations  were  completed,  the  decision  was  made  to  locate  the  FEL  axis 
57  inches  above  the  floor  Instead  of  58  Inches.  This  corresponds  to 
y  *  0.229  m,  a  change  whicji  is  small  enough  to  be  Ignored  here. 
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All  bending  magnets  produce  momentum  dispersion  In  the  bend  plane.  The 
extraction  line  must  be  made  non-dl sperslve,  l.e.,  achromatic,  so  that  the 
beam  can  be  properly  transported  through  the  FEL.  The  principle  of  the  method 
used  to  make  the  extraction  line  achromatic  Is  Illustrated  schematically  In 
figure  2.  In  figure  2,  we  show  a  beam  line  containing  two  dipole  magnets,  D12 
and  D13,  which  bend  through  equal  angles  In  opposite  directions  to  produce  a 
parallel -offset  beam.  The  dashed  line  Indicates  the  trajectory  of  a  particle 
whose  energy  Is  slightly  higher  than  the  central  trajectory  (which  Is  indica¬ 
ted  by  the  long-dash  —  short-dash  line).  The  high-energy  particle  would  be 
displaced  from  the  central  trajectory  by  an  amount  x  *  R16  (^)  after  leaving 
013.  Is  the  fractional  momentum  difference  between  the  high  energy  parti¬ 
cle  and  the  central  trajectory.  If  the  focusing  lenses  Q1  and  Q2  are  turned 
on,  the  high-energy  particle  will  follow  the  trajectory  Indicated  by  the  solid 
line  in  figure  2.  It  is  easy  to  see  that  R1S  will  be  zero  If  the  distances 

from  the  center  of  012  to  Q1  and  the  distance  from  Q2  to  the  center  of  D13  are 

S/4,  and  the  focal  lengths  of  Q1  and  Q2  are  S/8.  The  same  transformation  will 

make  R26  =  0,  where  x'  =  R26 

In  the  actual  situation,  we  have  dispersion  in  both  planes.  Since  all 
vertical  dispersion  originates  at  012V,  while  there  are  contributions  to  the 
horizontal  dispersions  from  several  magnets  (Dll,  El,  08,  and  D12H),  the  dis¬ 
persion  ray  Is  not  coplanar  with  the  system  axis.  This  is  not  really  a  pro¬ 
blem,  since  the  focusing  elements  are  quadrupoles.  Thus  each  thin  lens  in 
figure  2  Is  replaced  by  a  quadrupole  doublet  in  the  actual  system  of  figure  1. 
(Qx  ♦  Q12  and  Q13,  Q2  ♦  Q16  and  Q17).  The  system  can  be  made  achromatic  in 
both  planes  simultaneously.  This  solution  Is  described  In  table  II,  which 
gives  quadrupole  strengths,  beam  sizes,  and  acceptances  for  both  29  MeV  and 
185  MeV  beams. 
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The  only  important  shortcoming  of  this  "two  quad  pair"  solution  to  the 

achromatlclty  requirement  Is  that  beam  sizes,  particularly  In  Q16  -  Q17  and 

013  are  uncomfortably  large,  especially  at  low  extraction  energy.  At  29  MeV, 

the  system  acceptance  Is  l^ss  than  the  desired  value,  e  .  The  reason  for 

max 

this  is  easily  understood  from  figure  2.  The  two-dimensional  transfer  matrix 
between  D12  and  013  Is 


R 


two  quad 


r  i 


o' 

1 


(5) 


With  S  =■  8  m  determined  by  the  constraints  of  the  extraction  geometry, 

a  2  mr/mm,  meaning  that  an  Input  beam  size  of  1  mm  is  growing  at 

two  quad  21 

the  rate  of  2  mm  per  meter  of  path  length.  This  shortcoming  can  be  improved 
by  the  Insertion  of  a  third  quadrupole  doublet.  A  thin  lens  (QA)  Inserted  at 
the  center  of  the  schematic  system  of  figure  2  has  no  effect  on  the  dispersion 
but  changes  the  two  dimensional  transfer  matrix  to 


three  quad 


.  1 
S  f 


+1 


(6) 


where  f  is  the  focal  length  of  QA.  In  principle,  the  extraction  system  from 
D12  to  013  becomes  a  unit  matrix  transformation  when  f  *  S/16.  In  practice 
this  is  difficult  to  achieve  In  both  planes  simultaneously  and  requires  very 
strong  quadrupoles.  In  the  TRANSPORT  solution  to  the  "three  quad  pair" 
system,  we  used  the  central  quadrupole  doublet  Q14  -  Q15  to  minimize  the 
vertical  beam  size  at  013  after  observing  that  this  was  always  the  limiting 
aperture  of  the  system.  This  corresponds  to  f  *  S/8  In  the  schematic  model, 
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which  reduces  R21  by  a  factor  of  two  and  provides  adequate  aperture  size 
everywhere.  This  solution  is  described  In  table  III. 

We  will  implement  the  three-quad  pair  solution  to  be  certain  that  the 
system  has  the  required  acceptance.  The  two-quad  pair  solution  can  be  used  if 
the  design  goal  emittance  is  achieved.  It  will  be  possible  to  operate  the 
system  In  the  two-quad  pair  mode  by  turning  off  Q14  -  Q15  and  changing  the 
settings  of  the  other  quadrupoles,  since  the  locations  of  all  components  are 
the  same  for  both  solutions.  The  advantage  of  the  two  quad  pair  solution  is 
simplicity  -  fewer  adjustments  to  make  and  elimination  of  steering  corrections 
due  to  misalignment  of  Q14  -  15. 

IV.  PEL  Transport  Section 

The  FEL  transport  section  is  shown  schematically  in  figure  3.  Dipole 
magnets  D14,  15,  16,  and  17  comprise  a  chicane  whose  purpose  is  to  transport 
the  electron  beam  around  the  upstream  optical  cavity  mirror,6  indicated  by  the 
large  X  in  figure  3.  This  chicane  Is  manifestly  achromatic.  The  dump  magnet 
deflects  the  beam  through  90°  in  front  of  the  downstream  optical  cavity 
mirror.  Its  large  dispersion  provides  the  capability  of  energy  spectrum 
measurement.  All  bends  are  in  the  horizontal  plane.  The  two  quadrupole  pairs 
Q18-Q19  and  Q20-Q21  form  a  variable  magnification  telescope,  capable  of 
placing  beam  waists  of  adjustable  sizes  at  the  wiggler  center  in  both 
transverse  planes.  An  intermediate  waist  is  formed  by  the  telescope  in  the 
chicane  region,  which  tends  to  minimize  the  dipole  aperture  requirement. 

FEL  gain  is  maximized  by  maximizing  the  overlap  of  the  electron  beam  with 
the  optical  mode  in  the  cavity.  The  optical  mode  has  a  Gaussian  transverse 
intensity  profile  with  root  mean  square  radius  at  the  waist  given  by 
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*1 


(7) 


where  X  Is  the  optical  wavelength  and  z„  the  "Rayleigh  length."  The  root  mean 
square  optical  mode  size  Is  a  function  of  distance,  z,  from  the  waist,  given 
by 

o{z)  *  a0  /  1  +  (z/z0)z  .  (8) 


In  a  low-gain  FEL  system,  the  optical  waist  Is  placed  at  the  wlggler  center 
and  the  optimum  Rayleigh  range  Is  one  half  of  the  wlggler  length.  Me  maximize 
overlap  with  the  electron  beam  by  putting  Its  waist  at  the  wlggler  center, 
with  a  waist  size 


xo 


(9) 


where  L  is  the  wiggler  length  and  e  the  (unnormalized)  emlttance.  The 
electron  beam  size  at  the  ends  of  the  wlggler  will  then  be 

x  (1/2)  »/2x0  ,  (10) 

while  the  optical  mode  radius  is 

o(L/2)»/2a0  .  (11) 

Note  the  electron  beam  will  be  "Inside"  the  optical  mode  when 


e 


(12) 


The  situation  In  the  y-z  plane  Is  more  complicated  because  the  planar 

/ 

magnetic  wlggler  focusses  the  electrons  In  the  y-dlrectlon.7  (The  wiggle  is 

4-* 
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In  the  x  z  plane.)  The  y-dlrection  first  order  transfer  matrix  of  a  vrfggler 
segment  of  length  z  Is 


*  cos(j-) 

-  f-  sln(f-) 


sln(-) 

cos(J-) 


where 


"  2*  R 


is  the  wiggler  period  (28  mm),  and 

e  B0  X 

K  -  ~ z — “  (15) 

2/2  irmc 

is  the  wiggler  parameter  averaged  over  a  wiggler  period  (In  MKS  units).  This 
result  is  obtained  by  averaging  over  the  rapid  small  amplitude  motion  at  the 
wiggler  period.  In  the  same  approximation,  all  remaining  elements  of  the  six 
dimensional  transfer  matrix  are  given  by 


Rij  *  ■ 


except  R12  *  Z. 

In  our  TRANSPORT  calculations  we  represent  the  wiggler  by  four  drifts  of 
length  L/4,  plus  the  parameters  of  Ry  in  equation  (13),  supplied  numerically 
by  the  "Arbitrary  Transformation"  (type  code  14)  input  option.  We  require  a 
vertical  waist  at  the  wiggler  center  with  the  matched  beam  size 
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The  y  beam  size  then  remains  constant  within  the  wlggler,  as  can  be  calculated 
by  application  of  equation  (13). 

In  table  IV  we  present  a  summary  of  magnetic  field  strengths,  beam  sizes, 

and  ratios  of  planned  aperture  sizes  to  beam  sizes  for  a  185  MeV  electron 

beam.  Four  cases  are  presented:  (a)  two  quad  pair  extraction  system,  no 

wlggler  focusing;  (b)  same  Input  as  (a)  but  with  maximum  strength  wlggler 

focusing;  (c)  three  quad  pair  extraction  system,  no  wlggler  focusing;  and 

(d)  same  Input  as  (c)  but  with  wlggler  focusing.  Table  V  Is  constructed  like 

table  IV  for  a  29  MeV  beam.  Since  the  cases  based  on  the  two  quad  pair 

extraction  system  give  acceptances  much  less  than  e  „  ,  these  cases  will 

max 

probably  not  be  useful. 


Table  I  Extraction  Line  Layout 

r 

A.  Magnet  location  coordinates  [(a)  185  MeV,  (b)  29  MeV] 


Beam  Positions  in  Meters 

Magnet  Axis  Direction 

In  Degrees 

Entrance 

|  Exit 

- 

X 

y 

,  2 

y 

X 

z 

0x 

9y 

f 

Dll 

1.2359  (a) 
0.1730  (b) 

0 

0.950 

1.2300  (a) 
-0.1672  (b) 

0 

0.650 

180 

0 

0 

D8 

-0.0452  (a) 
-0.0658  (b) 

0 

0.275 

-0.0508  (a) 
-0.0695  (b) 

0 

0.425 

0 

0 

0 

D12H 

-0.0726  (a) 
-0.0762  (b) 

0 

1.050 

-0.0959 

0 

1.3494 

-2.000 

0 

0 

D12V 

-0.1079 

0 

1.4486 

-0.1440 

.0049 

1.7465 

-6.910 

0 

90 

Q12-Q13 

-0.3280 

.0551 

3.2645 

-0.3641 

.0650 

3.5622 

-6.910 

1.880 

0 

Q14-Q15 

-0.5806 

.1240 

5.3482 

-0.6166 

.1339 

5.6458 

-6.910 

1.880 

0 

Q16-Q17 

-0.8270 

.1913 

7.3823 

-0.8631 

.2011 

7.6798 

-6.910 

1.880 

0 

D13 

-1 .0388 

.2491 

9.1301 

-1  .0570 

.2540 

9.4299 

0 

0 

-15.260 

B.  Dipole  fields,  bend  angles  (a)  and  edge  angles  (01 ,  e2) 


Table  II  Extraction  Line  with  Two  Quad  Pairs 
A.  Quadrupole  Strengths 

Gradient  (gauss/cm) 


Quad  • 

29  MeV 

185  MeV 

Q10 

-  23 

0 

Qll 

22 

0 

Q12 

104 

637 

Q13 

-111 

-696 

Q16 

-113 

-705 

Q17 

116 

736 

Envelope  radii  (mm) 

29  MeV 

185  MeV 

Location 

x  y 

X 

y 

D12H-D12V 

1.81  2.24 

1.15 

1.11 

Q12-Q13 

1.77  5.10 

0.82 

2.00 

Q16-Q17 

10.27  8.42 

4.13 

3.80 

013 

2.86  2.90 

1.28 

1.39 

Limiting  Relative  Apertures* 

Aperture/beam  size 

29  MeV 

2.3 

185  MeV 

5.8 

Normalized  acceptance,  AN  (urn) 

26.5 

165 

♦Limit  occurs  in  Q16-Q17  vacuum  chamber,  assumed  physical  aperture  47.5  mm 
diameter 


/ 
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Table  III  Extraction  Line  with  Three  Quad  Pairs 
A.  Quadrupole  Strengths 

Gradient  (gauss/cm) 


Quad  ' 

29  MeV 

185  MeV 

Q10 

-  27 

0 

Qll 

27 

0 

Q12 

147 

871 

QI3 

-135 

-780 

Q14 

106 

763 

Q15 

-121 

-909 

Q16 

84 

533 

Q17 

-  81 

-588 

Envelope  radii  (mm) 

29  MeV 

185  MeV 

Location 

x  y 

X 

y 

012H-D12V 

2.23  1.98 

1.15 

l.n 

Q12-Q13 

1.15  5.89 

0.82 

2.32 

Q14-Q15 

4.19  5.90 

2.48 

1.53 

Q16-Q17 

3.28  1.97 

1.07 

0.42 

013 

1.51  2.41 

0.93 

1.08 

Limiting  Relative  Apertures 

Aperture  location 

29  MeV 

Q14-Q15(y) 

185  MeV 

D13(y) 

Assumed  aperture  size 

47.5  mm  dia 

20  mm 

Aperture/beam  size 

/ 

Normalized  acceptance. 

4.0 

9.3 

An  (ym) 

81.0 

430 

4.' 


Table  IV  FEL  Transport  Section  at  185  MeV 


Extraction  line 
focussing  type 

2  quad  pair 

2  quad  pair 

3  quad  pair  3  quad  pi 

Higgler 
f ocussl ng 

None, 

Maximum 

None 

Maximum 

A.  Quadrupole  Strengths  (gauss/cm) 

Q18 

-723 

-725 

-455 

-451 

Q19 

625 

625 

472 

469 

Q20 

915 

927 

637 

597 

Q21 

-526 

-549 

-591 

-490 

B.  Envelope  radii 

(mm) 

Location 

x  y 

x  y 

X 

y  x 

y 

Q18-Q19 

5.18  1.93 

5.19  1.93 

1.98 

1.20  1.98 

1.20 

Chicane 

3.35  1.08 

3.35  1.04 

1.18 

0.61  1.18 

0.67 

Q20-Q21 

0.39  0.26 

0.38  0.15 

0.47 

0.29  0.44 

0.20 

Higgler 

0.226  0.225 

0.224  0.143 

0.241 

0.240  0.240 

0.180 

Dump  Magnet 

0.27  0.34 

0.27  0.21 

0.30 

0.38  0.38 

0.26 

C.  Limiting  Relative  Apertures* 

Aperture/beam  size 

3.6 

3.6 

10.2 

10.1 

Normal  1  zed 

64.1 

64.2 

518 

514 

acceptance,  (urn) 

*L1m1t  occurs  In  014  vacuum  chamber.  In  the  bend  plane,  assumed  physical 
aperture  24.0  mm 
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Table  V  FEL  Transport  Section  at  29  MeV 


Extraction  line 
focussing  type 

2  quad  pair 

2  quad  pair 

3  quad  pair  3  quad  pal 

Wlggler 

focussing 

None* 

Maximum 

None 

Maximum 

A.  Quadrupole  Strengths  (gauss/cm) 

Q18 

-118 

-117 

-102 

-102 

Q19 

101 

101 

91 

91 

Q20 

174 

168 

109 

120 

Q21 

-138 

-167 

-108 

-147 

B.  Envelope  radii  (mm) 

Location 

x  y 

x  y 

X 

y  x 

y 

Q18-Q19 

12.33  4.05 

12.32  4.05 

3.95 

3.15  3.94 

3.15 

Chicane 

7.98  2.05 

7.89  2.15 

2.34 

1.63  2.27 

1.71 

Q20-Q21 

0.89  0.61 

1.13  0.41 

0.97 

0.54  1.15 

0.43 

Wlggler 

0.435  0.415 

0.398  0.143 

0.399 

0.398  0.398 

0.144 

Dump  Magnet 

1.08  1.66 

VO 
«— 1 

• 

CSJ 

o 

•-4 

• 

1.07 

1.36  1.10 

2.22 

C.  Limiting  Relative  Apertures* 

Aperture/beam  size  1.5  1.5  5.4  5.3 

Normalized  11.3  11.6  145  139 

acceptance,  AN  (um) 

♦Limit  occurs  in  014  vacuum  chamber,  in  the  bend  plane,  assumed  physical 
aperture  24.0  mm 
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APPENDIX  F 


Design  Note  5 

Predicted  FEL  Performance  at  74.375  MHz 


Predicted  PEL  Performance  at  74.375  MHz 
FEL  Design  Note  No.  5 

P.  Debenham,  C.  Johnson  and  S.  Penner 
July  20,  1987 


Hie  simplest  way  to  inject  electrons  into  the  RTM  for  FEL  operation  is  to 
illuminate  a  photocathode  directly  with  a  drive  laser.  This  will  provide  electron 
pulses  at  the  drive  laser  frequency,  74.375  MHz.1  Because  of  its  simplicity,  this 
mode  of  operation  is  especially  attractive  for  initial  operation  of  the  FEL.  The 
FEL  light  output  for  this  case  is  estimated  in  what  follows. 

The  small  signal  power  gain  of  the  NBS  FEL  has  been  calculated  by  Tang  et  al.2 
for  a  peak  electron  beam  current,  I,  of  2A.  It  is  8%  or  more  over  the  wavelength 
range  0.25 ^m  <  A.  <  10  ^m.  Assuming  a  power  reflectivity,  R,  of  99%  for  the  highly 
reflective  end  mirror  of  the  FEL  optical  cavity,  and  an  unspecified  additional  power 
loss,  u,  of  1%  per  cavity  round  trip,  we  can  choose  a  power  outcoupling,  C,  of  3%  to 
limit  the  round  trip  power  loss  to  6%.  This  leaves  a  2%  margin  to  the  gain  and 
ensures  that  the  FEL  will  lase. 

The  average  output  power  of  the  FEL  is  given  by  the  expression 


C 


Here,  N  is  the  number  of  periods  in  the  wiggler,  and  Pb  is  the  average  power  in  the 
electron  beam,  given  by 


Pb(W)  -  fi(MHz)I(A)r(ps)T(MeV), 


where  f^, z  and  T  are  the  pulse  frequency,  pulse  length  and  kinetic  energy,  respec¬ 
tively,  of  the  electron  beam.  T  may  be  written  as 

T  =  (^-DmQC2  =^mQC2. 

The  fundamental  wavelength, A  ,  is  given  by 


where  A w  is  the  wiggler  wavelength,  K  =  eB0Aw/27rm0c  (in  MRS  units),  and  Bq  is  the 


F  -  1 


peak  magnetic  field  in  the  wiggler.  Substituting  the  above  expressions  and  the 
following  values, 

C  *  0.03  £»  3  ps  N  ■  130  for  A<  2^m 

U  *  0.01  Aw  ■  2.8  cm  N  »  65  for  A  >  2 ji m 

R  =  0.99  B0  »  0.54  T  fj_  ■  74.375  MHZ 

I  -  2  A, 

we  obtain  a  formula  for  the  output  power  as  a  function  of  the  fundamental  wave¬ 
length: 

PoutW  *  73.3A"-5(^m)  for  \<  2pxn, 

*  147.  AT'5(/»*m)  for  A>  2^*m. 

As  is  shown  on  the  accompanying  graph,  this  mode  of  FEL  operation  can  provide 
between  45  W  and  165  W  of  radiation  over  the  full  wavelength  range  of  0.25 /u.m  to 

10  vi  m . 
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FEL  Design  Note  6 
FEL  Optical  Cavity 


C.  Johnson 
September  9,  1987 


I.  Optical  Resonator  Geometry 

The  inverse  of  the  round  trip  light  travel  time  of  the  optical  resonator 
is  constrained  to  be  a  submultiple  of  the  RTM  rf  frequency,  and,  in  addition, 
the  length  must  be  compatible  with  the  physical  dimensions  of  the  wiggler 
room.  The  design  value  is  L  «  806.158  cm.  The  Rayleigh  length,  zQ,  of  the 
optical  radiation  field  is  equal  to  half  the  length  of  the  wiggler,  and  the 
optical  beam  waist  is  centered  in  the  wiggler.  The  radii  of  curvature  of  the 
cavity  mirrors  are  constrained  by  L,  zQ,  and  the  position  of  the  optical 
waist.  The  parameters  for  the  optical  resonators  are  summarized  in  Table  1 
for  the  full  and  half  length  wigglers.  The  upstream  mirror,  M2,  is  highly- 
reflecting,  and  the  downstream  mirror,  Mj,  the  output  coupler,  is  partially 
transmitting.  Both  mirrors  are  formed  from  multi-layer  dielectric  coatings  on 
polished  substrates.  The  parameter  gi  =  1  -  L/R i ,  where  i  =  1,2  is  commonly 
encountered  in  optical  resonator  nomenclature;  the  product  gtg2  is  a  measure 
of  the  stability  of  the  cavity  (  jg^  |  <  1).  The  distance  t  is  from  the  cavi¬ 
ty  waist  to  Mj.  No*.e  that  the  optical  beam  radius,  oQ,  at  the  waist  is  given 
by  it 0^  =  \z  .  The  radius  scales  as  =  (1  +  1 4)1^. 

0  0  O  Z4, 


G  -  1 


4.’ 


Table  1 

NBS/NRL  FEL  Optical  Resonators 

Parameter  Full  Length  Wlggler  Half  Length. Wiggl  er 


Wavelength,  cm 

2.10-5  to  2.10-** 

2.10-4  to  lO"3 

Mirror  separation,  L,  cm 

Rayleigh  length,  zQ,  cm 

Cavity  waist  to  Mlt  t,  cm 

806.158 

182.0 

403.079 

806.158 

91.0 

494.079 

Outcoupler  radius,  Rlt  cm 

High  reflector  radius,  R2,  cm 

485.256 

485.256 

510.839 

338.614 

Cavity  stability  gt 

Parameters  g2 

9i92 

-0.661 

-0.661 

0.437 

-0.578 

-1.381 

0.798 

TEMoq  radius  oM  /O0 
on  mirrors  aJ,oQ 

(l/ez)  power  1 

2.430 

2.430 

5.521 

3.572 

Minimum  waist,  a0(^min)» 

0.340 

0.761 

Longitudinal  mode  spacing, 
c/2L  *  FSR,  MHz 

18.59 

18.59 

When  light  is  coupled  out  of  an  optical  cavity  through  a  partially  trans¬ 
mitting  coating  and  a  substrate  with  an  index  of  refraction  n(X),  the  sub¬ 
strate  acts  as  a  negative  lens,  so  that  the  Gaussian  beam  parameters  of  the 
emergent  TEMqo  mode  are  altered.  The  beam  waist  for  the  emergent  mode  is 
located  a  distance  h  behind  the  output  coupler  and  the  new  Rayleigh  length, 

z  is  smaller  than  z  —  the  emergent  beam  is  more  divergent.  The  waist 
o  o 

location  and  the  Rayleigh  length  depend  weakly  on  X  through  n(X).  The  sub¬ 
strates,  chosen  for  their  resistance  to  y  and  x-radiation,  are  UV  grade  fused 
silica  (S102)  (Suprasil)  for  0.2  urn  to  1.0  ym,  IR  grade  fused  silica  (Optosil  ) 
for  1.0  ym  to  2.5  ym,  and  zinc  sel  enlde  (ZnSe)  for  2.5  ym  to  10  ym.  Repre¬ 
sentative  values  for  h  and  z  1  are  given  in  Table  2. 

o 
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Table  2 


Emergent  TEMQ0  Mode 


X ,  urn 

n(X) 

h,  cm 

V’  cm 

Substrate 

0.2 

1.55051 

288.5 

84.0 

Suprasll 

0.3 

1.48719 

298.8 

90.7 

Suprasil 

0.4 

1.46962 

301.7 

92.7 

Suprasll 

0.5 

1.46252 

302.9 

93.5 

Suprasll 

0.6 

1.45840 

303.6 

94.0 

Suprasll 

0.7 

1.45515 

304.2 

94.4 

Suprasil 

0.8 

1.45330 

304.5 

94.6 

Suprasil 

0.9 

1.45170 

304.8 

94.8 

Suprasll 

1.0 

1.45024 

305.0 

95.0 

Opto si 1 

2.0 

1.43809 

307.2 

96.4 

Optosil 

2.75 

2.44 

208.2 

15.7 

ZnSe 

5.0 

2.43 

209.0 

15.8 

ZnSe 

7.5 

2.42 

209.9 

16.0 

ZnSe 

9.5 

2.41 

210.7 

16.1 

ZnSe 

11.0 

2.40 

211.6 

16.2 

ZnSe 

II.  Optical  Resonator  losses 

The  source  of  losses  considered  here  include  the  outeoupled  radiation, 
absorption  and  scattering  on  the  cavity  mirrors,  and  diffraction  losses.  I 
have  designed  the  cavity  such  that  the  apertures  of  the  wlggler  vacuum  chamber 
are  the  limiting  factor,  with  respect  to  diffraction  losses.  I  ignore  distri¬ 
buted  losses,  such  as  absorption  and  scattering  in  the  gain  medium.  The  out- 
coupled  fraction  is  variable,  T  =»  a  few  percent.  The  absorption  and  scatter¬ 
ing  losses  on  the  mirrors  are  taken  to  be  A  +  S  *  1  -  R^.  These  losses  are 
between  0.1%  and  5%  for  multi-layer  dielectric  mirrors.  (Typical  values  for 
R^r  may  be  found  from  manufacturer's  data.) 

Diffraction  losses  were  estimated  from  published  solutions1  to  the  inte¬ 
gral  equations  that  result  from  the  Fresnel  -Kirchhoff  formulation  of  Huygen's 


lH.  Kogel  nik  and  T.  LI,  Proc.  IEEE,  54,  1312  (1966),  and  G.D.  Boyd  and  J.P. 
Gordon,  Bell  Sys.  Tech.  0.,  40,  489  TT961). 
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principle.  For  the  case  of  square  mirrors  (or  apertures),  the  field  distribu¬ 
tions  are  given  approximately  by  Hermite-Gaussian  functions,  and  for  the  case 
of  circular  mirrors,  approximately  by  Laguerre-Gaussian  functions..  The  losses- 
were  calculated  only  for  (the  TEMQ0  mode.  In  the  case  of  an  asymmetric  aper¬ 
ture,  the  smaller  dimension  was  taken,  so  these  diffraction  losses  are  conser¬ 
vative  estimates.  Apertures  were  treated  as  if  they  were  mirrors  with  a  radi¬ 
us  of  curvature  equal  to  that  of  the  wavefronts  of  the  optical  field  at  the 
location  of  the  aperture. 

Kogelnik  and  Li  publish  solutions  for  symmetric  (gj  =  g2)  cavities  with 
circular  mirrors  for  the  confocal  case  (|g|  *  0)  through  to  the  concentric/ 
planar  case  (|g|  =  1).  Boyd  and  Gordon  publish  solutions  for  square  mirrors 
and  symmetric,  confocal  geometry.  We  have  a  symmetric  situation  for  the  full 
length  wiggler  end  mirrors  and  wiggler  vacuum  chamber  apertures,  but  otherwise 
the  situation  is  asymmetric.  In  all  cases  the  geometry  is  non-conf ocal . 

Asymmetric  resonators  were  transformed  into  symmetric  resonators  with  the 
same  diffraction  loss  using  the  equivalence  principles  for  spherical  mirror 
resonators: 

1.  The  mirrors  can  be  interchanged; 

2.  The  sign  of  both  of  the  stability  parameters  can  be  reversed,  gt,  g2 
=  ”9 i »  “92* 

3.  Resonators  with  the  same  values  for  N,  Gj,  and  G2  are  equivalent, 
where  N  is  the  Fresnel  number,  N  *  a1a2/(\L),  Gt  =  9iai/a2»  and 
G2  =  g2a2/ai* 

Here  2 a  is  the  diameter  of  a  circular  mirror  or  the  length  of  a  square  mirror. 

In  other  words,  using  principle  3,  a  resonator  with  g1  *  g2  was  transformed 
into  a  resonator  with  Gj  =  G2  by  varying  the  ratio  of  at  to  a2.  Then,  for 
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circular  mirrors,  the  diffraction  loss  per  reflection  was  calculated  by  using 
the  value  of  S  for  g,  calculating  N,  and  using  Figure  22  In  Kogelnlk  and  Li, 
For  an  asymmetric,  non-confocal  cavity  with  square  mirrors,  I  estimated 
the  diffraction  losses  by  transforming  to  a  symmetric  cavity  as  before  and 
then  assuming  that  the  loss  Is  equal  to  that  of  a  confocal  resonator  with  the 
reflector  dimensions  scaled  by  the  ratio  of  the  spot  sizes.  This  amounts  to 
Mderat1ng“  the  Fresnel  number  by  the  factor 
N*  *  N(2L/R  -  (L/R)2)1/2, 
and  then  using  Figure  3  In  Boyd  and  Gordon, 

The  results  of  the  calculations  are  given  In  Table  3,  which  lists  the 
diffraction  loss  per  reflection,  Dv,  for  the  significant  apertures.  The  point 
z  *  0  corresponds  to  the  waist  of  the  resonator.  Note  that  for  efficient 
operation  of  the  half  length  wiggler,  the  vacuum  aperture  for  the  dump  magnet 
should  be  Increased  to  2.5  x  2.5  cm. 


III.  FEL  Operating  Characteristics 

In  order  to  lase,  the  losses  must  be  smaller  than  the  small  signal  power 
gain.  By  estimating  the  losses  and  the  small  signal  power  gain.  I've  deter¬ 
mined  the  lower  limit  for  the  wavelength. 

I  calculated  the  small  signal  power  gain  from 


G  3  f2  iL  J_  |t2m3  i_  sin^v 

p  1  !A  ^  3v  ^  ’ 

3  sin2v  r, 

3V^ - *  0.54, 


3  s1n2v 
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I 

•> 

8  *  5400  Gauss,  so  that  K  =  1.41,  and  then  a 

N  3  130  or  65  (full  length  or  half  length  wiggler), 

Atl  =  2.8  cm 
w 

1=2  Amps,  and 

I 

tro*  =  AZ  . 

0  0 

By  substituting  for  y2  from  the  FEL  resonance  relation. 


I  get  Gp(A)  3  17.09  /A  full  length  wiggler 
Gp( x)  =  4.272  /A  half  length  wiggler, 
where  A  is  in  cm.  This  formula  for  the  small  signal  power  gain  is  based  on 
the  1-D  calculations  and  with  a  conservative  estimate  for  the  filling  factor 
(the  overlap  between  the  electron  and  optical  spatial  modes  inside  the  wig¬ 
gler).  C.M.  Tang  has  revised  the  estimation  for  the  small  signal  gain  to 
include  several  3-0  effects,  and  a  more  realistic  estimation  for  the  filling 

factor.  The  overall  result  is  about  the  same  values  for  G„  as  what  I  have 

P 

used. 

According  to  the  engineers  at  Acton  Research,  it  is  not  possible  to  con¬ 
struct  multilayer  dielectric  coatings  for  wavelengths  below  146  nm.  Typical 
reflectivities,  taken  from  the  Acton  Research  catalog,  at  146  nm  are  RHR  = 
0.94  t  0.02.  At  2  Amps  peak  current,  we'll  have  Gp ( 146  nm)  s  6.5%,  which 
requires  R^R  >  0.97  (see  below).  Below  146  nm,  metallic  coatings  are  used, 
with  typical  reflectivities  of  0.70  to  0.90.  Therefore,  without  an  improve¬ 
ment  in  the  small  signal  power  gain  (or  better  mirrors),  lasing  below  146  nm 
will  not  be  possible  with  the  NBS/FEL.  Note  that  below  200  nm,  we  would  have 
to  change  to  MgF2  or  Li F  for  all  of  the  transmissive  optics.  Below  110  nm. 
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no  material  transmits,  and  this  defines  the  start  of  the  “extreme 
ultraviolet,"  or  EUV,  region. 


After  lasing  begins,  the  PEL  will  saturate.  The  output  no  longer  depends 

on  the  small  signal  power  gain.  In  theory,  the  average  power  increase  per 

\ 

round  trip  is  ^b/2N,  where  T5^  Is  the  average  power  in  the  electron  beam.  The 
beam  power  is  determined  by  the  kinetic  energy  of  the  electrons,  the  peak 
current,  the  pulse  width,  and  the  frequency: 

T^(a)  3  fj»I»At  (ymoc2)  Watts,  or 

( x ) / 2N  =  1.465//Y”  Watts  (full  length  wiggler) 

=  2.931  //T"  Watts  (half  length  wiggler), 
where  x  is  in  cm  and  I've  used  the  PEL  resonance  relationship  to  eliminate  y. 

Por  the  saturated  PEL,  it  is  straightforward  to  calculate  the  average 
intracavity  and  output  power  in  terms  of  F^fx).  After  one  complete  round 
trip,  the  power  in  the  optical  field  must  be  the  same,  i.e.,  the  fractional 
increase  in  average  power,  7k/(2NTT  ),  is  equal  to  the  sum  of  the  fractional 
losses  --  here  diffraction  at  the  vacuum  apertures  (D^ )  and  cavity  mirrors 
(D^),  absorption  and  scattering  losses  on  the  mirrors,  and  transmission  loss 
(the  output  power).  The  average  normalized,  intracavity  power  on  the  output 
cou pi er  is 


cav.oc 

Pb/2N 

The  average. 


11  -  0,) 

2  ‘  '  i  -  roc»  -  V  <l  -  V 

normalized  output  power  is 


yout  (1  -  °„H1  -  °,)T _ 

Pb/2N  1  ‘  RHR  R0C(1  ’  °v)H  ^  ‘  V 
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At  the  beginning  of  the  wiggler,  the  average,  normalized  intracavity  power  is 

‘  r%7-  "Iff  ■  "ocV1  -  0,i'»  -  V2- 

b' 

The  condition  that  the  small  signal  fractional  power  gain  must  exceed  the 
fractional  losses  can  now  be  made  quantitative: 

G  >  1  ~  Reff  =  _L_ 
p  *eff  em1n 

The  higher  RHR  is,  the  easier  it  will  be  to  achieve  saturation,  the  larger  T 
can  be,  and  the  more  efficient  the  PEL  will  be.  For  higher  values  of  Gp,  the 
minimum  value  of  RHR  decreases,  i.e.,  R^R  >  0.90  at  2  pm  for  the  full  length 
wiggler,  where  6p  =  24%.  If  absorption  and  scattering  losses  are  small  at 
this  wavelength,  say,  0.1%,  then  T  could  be  as  large  as  19%. 

Efficient  laser  operation  requires  low  losses.  However,  if  the  losses 
(including  T)  are  made  small  enough,  the  intracavity  power  will  reach  levels 
that  are  capable  of  mirror  damage.  Therefore,  we  must  consider  damage  thres¬ 
holds  of  the  multilayer  dielectric  coatings  when  chosing  the  appropriate  val¬ 
ues  for  R^r  and  T.  The  average  irrad lance  (W  cm-2)  and  peak  fl  uence  (J  cm-2) 
are  the  important  quantities,  so  we  must  calculate  the  effective  area  of  the 

radiation  on  the  cavity  mirrors.  The  area  is  Tra£  *  *a2  (— )2,  where  —  =  2.43 

M  o  o'  a 

o  o 

(full  length  wiggler)  and  —  =  3.57  (half  length  wiggler)*  (see  Table  1).  The 

°o 

average  intracavity  irradiance  and  peak  intracavity  fl  uence  on  the  output 
coupler  can  be  put  in  terms  of  the  PEL  wavelength  using  the  expression  for 
Fb(A)/2N: 


♦This  is  a  conservative  estimate,  because  I've  used  o  /cr  for  the  high  re¬ 
nt  o 


fl  ector. 
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gPK  1 

T  *  nr-  -  71-  (— )2,  W  cm"2 
oc  2N  Xz„  '  au'  * 

0  M 

F  =»  4^-  J  cm-2, 
oc  ri 


Numerically, 


I 


00  734«X3/2 


^575  W  cm"2  >  -*?8  •710"5  w  cm-2. 


e  emin 


c  -  e  •  1.83  »  10~n  .  2  i 
F  —  J  cm  ^ 

oc  ,3/2  ♦ 

X  e  =  e 


1.07  •  10-12  , 

- -2 - j  cm  S 


for  the  full  length  wiggler,  and 


I  = 


396*X 


3/2 


W  cm"2  > 


mi  n 


5.915  •  10"1* 


W  cm-2. 


£  ~  Emin 


c  _  e  •  3.39  •  10"11  ,  „  7.94  •  10‘12  , 

Foc - 372 - J  cm  l  - 77 - Jem2, 

X  e  =  emin 


T 


for  the  half  length  wiggler.  As  before,  the  wavelength  is  in  cm.  The  minimum 
intracavity  average  irradiance  on  the  output  coupler  is  plotted  in  Figure  1 
(solid  line),  along  with  observed  cw  damage  thresholds  for  high  power,  multi¬ 
layer  dielectric,  laser  mirrors  (dashed  line),  overcoated  silver  (dash-dot 
line),  and  overcoated  aluminum  (dotted  line).  If  we  regard  the  FEL  as  a  cw 
laser,  then  these  metallic  coatings  are  too  fragile  for  our  use,  because  they 
are  destroyed  before  the  1  aser  saturates.  The  individual  points  represent 
calculated  solutions. 

Damage  thresholds  for  pulsed  laser  radiation  are  generally  quoted  for  10 
ns  wide  pulses  at  10  Hz  repetition  rate  at  1.06  um.  To  derate  for  our  short 
pulse,  I  fitted  observed  breakdown  fluences  (see  Melles  Griot  Optics  Guide  3) 

^  n 

at  several  pulse  widths  to  the  form  FgD  «  t  to  find  FgQ  =  1.6  •  106  t  J 
cm"2,  or  Fgp  =  0.6  J  cm"2  for  3  ps  wide  pulses.  The  observed  data  were  for 
near-IR  wavelengths;  standard  practice  is  to  derate  by  a  factor  of  5  for  the 
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UV  and  by  a  factor  of  2  for  the  visible.  The  minimum  intracavity  fluence  on 
the  output  coupler  is  plotted  in  Figure  2,  along  with  the  estimated  damage 
thresholds  for  high  power  laser  mirrors  with  3  ps  wide  pul  ses,  and  the  calcu¬ 
lated  solutions. 

♦ 

With  an  idea  of  permissable  build-up  ratios,  typical  reflectivities,  and 
diffraction  losses,  we  can  estimate  the  output  power.  Note  that  to  convert  to 
output  irradiance  or  fluence,  the  correct  area  is  Az^  (see  Table  2).  The 
results  are  summarized  in  Table  4  and  in  Figures  3,  4,  and  5.  In  all  cases, 
the  pulse  width  is  3  ps  and  the  repetition  rate  is  74.375  MHz.  The  average 
diffraction  losses  are  given  by  U,  where  D  =  1  -  (1  -  0 v ) 2 ( 1  -  D^).  The 
efficiency  is  defined  for  the  overall  process,  taken  to  be  the  product  of  the 
RTM,  wiggler,  and  optical  cavity  efficiencies:  Efficiency  =  — — .  The 

parameter  e1/2  is  the  half-angle  of  the  divergence  of  the  optical  TEMqq  mode: 

ei/2  =  tan  °o/Zo  a  ao/zo  =  X/nao- 

In  addition  to  the  output  power  and  spatial  mode  characteristics,  there 
are  a  number  of  other  important  parameters.  The  output  is  linearly  polarized. 
The  spectral  bandwidth  is  expected  to  be  Fourier-transform  limited,  so  that 
the  bandwidth  depends  on  the  pulse  width  and  shape:  for  pulses  with  a  Gaus¬ 
sian  temporal  distribution  with  a  FWHM  of  3  ps,  the  spectral  bandwidth  is 
about  200  GHz,  or  7  cm-1,  or  1  .7  A  at  500  nm.  This  bandwidth  corresponds  to 
about  104  longitudinal  modes  of  the  8.06  m-1  ong  oscillator.  If  the  pulse 
width  is  3  ps  over  the  entire  wavelength  range,  then  the  resolution,  AX/A,  is 
the  smallest  in  the  UV.  It  should  be  possible  to  compress  the  pulse  width  to 
about  100  fs. 

The  tuning  range  of  standard  high-power,  multi-layer  dielectric,  laser 
mirrors  is  about  ±  12%;  so-called  broad  band  mul  ti-1  ayer  diel  ectric  reflectors 
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are  about  a  factor  of  two  broader  —  this  corresponds  to  a  bandwidth  of  200- 
250  n m  in  the  visible.  Hence,  in  order  to  oscillate  in  the  fundamental  from 

0.2  um  to  10  pm,  many  mirrors  will  be  needed.  However,  it  is  important  to 

realize  that  the  PEL  will  be  flexible  without  changing  mirrors  —  the  FEL 

I 

output  can  be  used  to  generate  harmonics  from  second -harmonic  generation  and / 
or  mixing  in  nonl  inear  crystals,  four  wave  mixing  in  gases,  and  Raman  shifting 
in  H2.  In  addition,  the  FEL  will  have  sufficient  power  to  pump  a  number  of 

tunable  lasers,  i.e.  organic  dye  lasers  for  the  near-UV  to  near  IR,  and  color- 

center  or  FIR  lasers  in  the  infrared.  It  may  be  possible  to  synchronously 
pump  a  dye  laser  in  order  to  generate  femtosecond  pulses. 

IV.  Cavity  Modifications 

There  are  several  disadvantages  to  the  use  of  partially  transmitting/ 
partially  reflecting  multi-layer  dielectric  coatings  for  the  output  coupler. 
The  bandwidth  may  be  narrower  than  that  of  the  high  reflector,  T  is  variable 
by  as  much  as  50%  for  some  coatings,  and  to  change  T  one  must  change  the 
mirror.  A  possible  alternative  to  this  cavity  would  consist  of  two  highly 
reflecting  end  mirrors  and  an  intracavity,  variable  output  coupler,  for  exam¬ 
ple  an  electrooptic  modulator  and  a  polarization  beamsplitter.  If  the  E/0 
crystal  were  cut  to  act  like  a  half-wave  plate  at  zero  applied  voltage,  then 
the  output  would  vary  as  sin2r,  where  r  is  small  and  proportional  to  the 
applied  voltage.  The  disadvantage  to  this  scheme  is  that  the  E/0  crystal  and 
polarization  beamsplitter  must  be  anti-refl  ection  (AR)  coated  to  minimize 
losses,  and  these  coatings  also  have  limited  bandwidth.  If  the  voltage  was 
switched  rapidly,  so  that  r  changed  from  0  to  ±  n/2,  then  the  E/0  modulator 
could  be  used  as  a  cavity  dumper. 
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An  alternative  method  of  cavity  dumping  involves  acousto-optic  modula¬ 
tion.  Basically,  a  pulse  of  rf  voltage  is  applied  to  a  transducer  that  is 
mounted  on  a  piece  of  optical  quality  quartz.  Acoustic  waves  are  generated, 
and  it  Is  possible  to  deflect  an  optical  beam  off  the  diffraction  grating  set 

t 

up  by  the  acoustic  waves.  The  A/0  cavity-dumper  is  superior  to  the  E/0  modu¬ 
lator  because  the  quartz  can  be  cut  at  Brewster's  angle,  thus  eliminating  the 
need  for  AR  coatings,  and  therefore  extending  the  tuning  range.  For  efficient 
use  of  the  A/0  beam  deflector,  it  is  necessary  to  have  an  optical  waist  in  the 
quartz,  so  that  a  three  mirror  FEL  optical  cavity  might  be  necessary. 
Acousto-optic  cavity  dumpers  can  be  used  as  output  couplers  in  conventional 
lasers.  An  application  involving  a  dye  laser  is  shown  in  Figure  6. 
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FIGURE  I 


AVERAGE  INTRACAVITY  1RRA0IANCE  ON  THE  OUTPUT  COUPLER 
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FIGURE  2 


PEAK  INTRACAVITY  FLUB'CE  ON  THE  OLTIPUT  COUP^R 
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SYNCHRONOUSLY  -  PUMPED.  CAVITY  -  DUMPED.  DYE  LASER:  SEE 
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FIGURE  6 


FEL  Design  Note  7 
C.  Johnson 


November  18,  1987. 

Tuning  the  optical  wavelength  and  Diffraction  losses  in  the  IR 

We  expect  to  vary  the  FEL  optical  wavelength  by  changing  the 
electron  beam  energy  and  the  value  of  the  peak  magnetic  field  of  the 
wiggler.  We  will  scan  the  optical  wavelength  while  maintaining 
approximately  constant  optical  power  by  varying  the  magnetic  field.  The 
electron  beam  energy  will  be  given  by  the  number  of  orbits  through  the 
accelerator  and  by  the  energy  gained  per  orbit.  For  a  ±  20%  variation 
per  orbit,  the  energy,  in  MeV,  is  given  by 

T  -  5  +  Norblt  (  12  ±  2.4  )  MeV. 

The  values  of  the  beam  kinetic  energy  determine  the  optical 
wavelength  and  the  small-signal,  one -dimensional,  fractional  power  gain 
Gp  in  the  usual  way  (.see  earlier  design  notes).  The  fiducial  parameters 
are  defined  in  Design  Note  6.  It  is  interesting  to  note  that  for  Norbit 
>  2,  the  same  optical  wavelength  can  be  obtained  for  two  or  more  values 
of  Norblt  (see  Figure  1).  In  fact,  at  B0  =  5400  Gauss  the  only  gap  in 
the  optical  wavelength  occurs  between  12  /jm  and  18.4  pm,  corresponding  to 
the  low  energy  limit  of  orbit  2  and  the  high  energy  limit  of  orbit  1 .  In 
theory,  it  is  possible  to  cover  this  gap  by  decreasing  the  peak  magnetic 
field  from  B0  to  0.6  B0  with  T  held  fixed  at  19.4  MeV  (e.g.,  see  Figure  2 
in  Tang  et.  al.,  PAC  Conf.  Proc.  1987)  ,  so  that  the  FEL  could  be  tuned 
from  150  nm  to  32  pm.  However,  diffraction  losses  in  the  IR  will  not 
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support  such  a  broad  dynamic  range  of  operation  (see  below) . 

To  estimate  diffraction  losses,  I  have  followed  the  procedure 
outlined  in  Design  Note  6  (see  D.N.  6,  Table  3),  but  with  the  wavelength 
extended  to  35  nm.  The  fractional  power  loss  per  aperture  is  plotted  in 
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Figure  2  for  the  wiggler  vacuum  chamber,  the  dump  and  D17  magnet  vacuum 
chambers ,  and  the  cavity  end  mirrors .  The  dump  magnet  vacuum  chamber  is 
taken  to  be  2.5  cm  x  2.5  cm.  The  quad  magnet  vacuum  chambers  do  not 
contribute  to  the  diffraction  losses. 

To  compare  these  losses  to  the  small-signal  one -dimensional  power 
gain,  I  assume  that  diffraction  is  the  only  loss  mechanism: 

R2di££  »  (  1-  Dmr  )2  (  1  -  D*  )*  (  1  -  dmag  )* 

where  MR ,  w,  and  MAG  stand  for  the  cavity  end  mirrors,  the  wiggler  vacuum 
chamber,  and  the  dump/D17  magnet  pair.  The  quantity  Rzdlff  is  compared 
to  the  minimum  value  that  is  required  to  overcome  the  oscillator 
threshold: 


mix,  1  +  Gp 

The  maximum  wavelength  is  around  12  ^m  (see  Figure  3). 

In  an  actual  optical  resonator,  the  reflectivity  of  the  cavity 
mirrors  would  have  to  be  included: 

R2«ff  =rhr  Roc  R2 d i t f  > 

so  that  the  maximum  value  of  the  optical  wavelength  would  be  less  than  12 
/Lim.  It  should  be  emphasized  that  the  diffraction  losses  in  this  design 
note  are  estimated,  and  that  the  small  signal  gain  is  for  the  one 
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dimensional  case.  It  may  be  that  optical  guiding  will  offset  diffraction 
so  that  the  FEL  will  lase  at  wavelengths  longer  that  12  nm. 

Erratum  to  Design  Note  6 

The  stability  condition  on  page  1  should  read  0  <  gt  g2  <  1  , 
not  |  g!  g2  |  <  1. 
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Design  Note  8 

Electron  Beam  Phase  Stabilization 


FEL  Design  Note  #8 
Electron  Beam  Phase  Stabilization 

Samuel  Penner 
November  4,  1987 

t 

In  a  recent  conversation  with  Tom  Baer,  he  asked  If  It  would  be  possible 
to  stabilize  the  electron  beam  phase  beyond  the  ±0.5°  that  we  expect  to 
achieve*  His  reason  for  doing  so  Is  that  If  the  beam  phase  and  therefore  the 
RF  frequency  is  sufficiently  stable.  It  might  (when  combined  with  several 
other  features)  allow  frequency  locking  between  the  RF  and  the  FEL  optical 
output.  This  could  have  very  Important  metrological  applications.  The  sta¬ 
bilization  process  would  require  a  large  bandwidth  because  some  very  high 
harmonic  of  the  RF  frequency  Is  needed  to  begin  the  chain  of  frequency  multi¬ 
plications  which  relates  RF  frequency  to  optical  frequency. 

A  preliminary  look  at  this  problem  Indicates  that  beam  phase  stabiliza¬ 
tion  Is  possible  by  a  feedback  scheme,  provided  that  a  phase  error  signal  can 
be  derived  from  the  beam.  The  feedback  device  is  an  RF  cavity  which  modulates 
the  beam  energy  preceding  a  portion  of  the  beam  transport  system  which  pro¬ 
vides  time  dispersion  (i.e.  transport  matrix  element  *56  *  0).  Specifically, 
a  cavity  with  peak  RF  voltage  AV  located  on  the  last  return  leg  of  the  RTM 
will  produce  a  phase  change  at  the  FEL  of  A+  ■  RS6(aV/E0)  where  E0  Is  the 
extracted  beam  energy.  For  the  case  of  extraction  at  185  MeV,  our  earlier 
transport  calculations  give  R56  *  -54  degrees  per  percent  energy  change.  If 
the  feedback  cavity  Is  phased  so  that  the  beam  traverses  It  at  Its  peak  phase, 
a  cavity  voltage  amplitude  of  34  kV  Is  required  to  change  the  beam  phase  by  1° 
at  the  FEL.  A  second  cavity  located  on  the  extraction  line  (s«y  on  the  seg¬ 
ment  between  Q19  and  D14)  with  the  same  amplitude  and  opposite  phase  relative 
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to  the  beam,  will  eliminate  the  energy  modulation.  The  effect  of  these  cavi¬ 
ties  on  the  transverse  optics  Is  negligible.  The  effect  on  the  energy  and 
phase  spread  of  the  beam  is  likewise  negligible  because  of  the  small  phase 
spread  In  the  beam.  In  principle,  by  operating  at  1  arger  amplitudes  and  dif¬ 
ferent  phases,  these  same  cavities  could  be  used  to  manipulate  the  beam  longi¬ 
tudinal  phase  space,  e.g.,  to  shorten  the  pulse  from  Its  nominal  3°  phase 
spread. 

A  single  2380  MHz  cavity  can  provide  the  34  kV  amplitude  drive.  It  must 
be  a  low  Q  cavity  because  we  would  like  to  have  feedback  response  at  the  FEL 
optical  cavity  frequency,  which  is  the  128th  subharmonic  of  the  RF  fundamental 
frequency  (18.594  MHz).  This  requires  a  cavity  Q  of  order  100.  Since  (R/Q)  * 
350  0  for  the  cavities  In  our  main  llnac  sections,  a  cavity  of  the  same  geome¬ 
try  with  Q  »  100  requires  33  kW  of  drive  power.  This  Is  high,  but  not  in¬ 
tractable.  The  power  required  could  be  further  reduced  by  using  a  multiple 
cell  cavity,  e.g.,  a  three  cell  cavity  would  require  only  11  kW.  It  may  also 
be  possible  to  use  less  cavity  vol  tage.  At  present  the  observed  beam  phase 
jitter  at  5  MeV  Is  less  than  1°  (peak  to  peak).  Some  of  this  jitter  origi¬ 
nates  in  noise  in  the  high  voltage  terminal  and  very  likely  can  be  reduced. 
Further,  if  the  phase  jitter  has  no  frequency  components  in  the  18  MHz  range, 
the  cavity  Q  can  be  Increased  to  reduce  the  power  needed. 

To  preserve  maximum  bandwidth,  the  sensor  which  generates  the  phase  error 
signal  on  the  beam  must  be  located  near  the  RTM.  It  would  not  be  feasible  to 
generate  the  phase  error  signal  at  the  FEl  because  the  transit  time  of  the 
beam  from  the  feedback  cavity  to  the  FEL  Is  approximately  80  ns.  It  would 
take  a  similar  time  for  an  error  signal  to  travel  back  to  the  drive  cavity. 
From  the  propagation  delays  alone,  the  maximum  feedback  bandwidth  must  be  less 
than  6  MHz. 
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While  It  would  be  Interesting  to  pursue  these  possibilities  further,  thi 
development  is  beyond  the  scope  of  the  present  FEL  project,  and  we  have  no 
funding  at  present  to  implement  it. 
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Lasing  on  the  Harmonics 


FEL  Design  Note  9 
Lasing  on  the  Harmonics 

C.  Johnson 


March  15,  1988 

The  short  wavelength,  or  ultraviolet,  cuto££  for  lasing  with  the  NBS/FEL 
is  determined  by  the  point  where  the  small  signal  power  gain  no  longer  exceeds 
the  cavity  losses.  This  design  note  shows  that  we  can  operate  at  shorter 
wavelengths  while  sacrificing  output  power  by  adjusting  the  electron  energy 
for  operation  at  At  ,  the  "fundamental  wavelength”,  and  adjusting  the  cavity 
mirrors  for  operation  at  A3  =  At/3,  the  "third  harmonic  wavelength".  This  is 
because  for  fixed  lasing  wavelength,  in  the  one-dimensional  approximation,  and 
with  the  maximum  magnetic  field,  the  small  signal  power  gain  in  the  third 
harmonic  is  higher  than  that  of  the  fundamental.  This  is  also  true  for  the 
fifth  harmonic,  but  not  for  harmonics  with  n  i  7. 

The  maximum  value  of  the  one - dimens ional ,  small-signal  power  gain,  Gn , 
for  the  nth  harmonic  is 


2  2  3 

G  =  n  F  (K)  K  N 
n  n 


t  V 
I  X  A 

W 


h  V 


3  ’ 


where  the  symbols  have  their  usual  meaning  [Ref.  1).  The  ratio  of  the  gain  in 
the  harmonic  to  that  of  the  fundamental  is 


Gn(K) 

G^K) 


1  F»(K)  ^ 
1  F1<K>  ^n 


f2(K) 
n  fJ(K) 


5/2 

n 


Fn<*> 

fJ(K> 


Here  A  is  the  lasing  wavelength.  The  cross  sectional  area  of  the  radiation 
field  cancels  since  the  lasing  wavelength  is  the  same  in  either  case.  The 
functions  Fn  involve  the  differences  of  Bessel  functions,  where  the  argument 
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depends  on  Che  harmonic  index  and  Che  value  of  Che  wiggler  paramecer,  K.  For 
Che  NBS  FEL,  K  is  beCveen  1.0  and  0.6  (rms).  SubsCiCuCing  for  Fn  [Ref.  2], 


Gn(K)  , 

g]7k)  * 


n5^ 


Jn-l<n<>  -  -Wn«>  I' 

T  ~ 


[Jo«>  *  Ji<?)]2 


where  £  =  Kz/(2(1+K2)) .  Ac  K  =  1.0,  i  =  0.25  and  G3/G1  =  1.677  and  Gj^  = 
1.360.  AC  K  =  0.6,  £  =  0.13  and  Gn  <  Gt  for  all  values  of  n. 

Having  esCablished  ChaC  Che  harmonic  gain  can  exceed  ChaC  of  Che  funda- 
mencal  in  cerCain  limiCs,  che  acCual  value  of  Che  gain  muse  be  deCermined. 

The  one-dimensional  formulaeion  neglecCs  several  effeecs  and  cannoC  be  used  Co 
predice  absoluCe  values  for  Che  small  signal  power  gain.  Assuming  ChaC 
Che  raCio  of  Che  gain  in  a  harmonic  Co  ChaC  of  Che  fundamencal  is  given  by  Che 
one -dimensional  resulc,  I  have  calculaCed  Che  harmonic  gain  from  Tang’s  3-D 
calculacions ,  wich  normalized  beam  emiccance  of  10  mm-mrad.  I  ficced  three 
points  in  her  figure  (see  Figure  6  in  Ref.  1),  corresponding  to  K  =  1.0,  en  = 
10  mm-mrad,  and  7  =  350,  250,  and  150,  to  the  function  alb  .  The  result  is 
G1(3-D)  =  31.71°  71 ,  where  X  is  in  pm.  Then,  the  gain  in  the  third  harmonic 
is  1.68  times  G1(3-D).  This  information  is  plotted  in  Figure  1,  along  with 
representative  values  for  the  overall  cavity  power  loss  per  pass.  The  cavity 
losses  were  calculaCed  using  excimer  laser  mirror  specif ications ,  see  Table  1. 
The  outcoupling  factor  is  5%  for  all  wavelengths.  The  ability  to  lase  on  the 
third  harmonic  extends  our  operating  range  from  200  nm  to  about  150  nm. 
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Table  1 


Wavelength 

Mirror  Parameters 

Power  Reflection  Coefficient 

i— t-B2.: 

nm 

High  Reflector 

Output  Coupler 

p2»f  f 

146 

0.960 

0.910 

0.145 

158 

’  0.970 

0.920 

0.121 

172 

0.980 

0.930 

0.097 

193 

0.990 

0.940 

0.075 

212 

0.990 

0.940 

0.075 

222 

0.990 

0.940 

0.075 

248 

0.990 

0.940 

0.075 

266 

0.990 

0.940 

0.075 

282 

0.990 

0.940 

0.075 

308 

0.990 

0.940 

0.075 

325 

0.999 

0.949 

0.055 

337 

0.999 

0.949 

0.055 

352 

0.999 

0.949 

0.055 

400 

0.999 

0.949 

0.055 

1000 

0.999 

0.949 

0.055 

In  the  low-gain,  one -dimensional  approximation,  the  coupling  efficiency 
rj ,  which,  in  the  case  of  an  untapered  wiggler  and  a  lossless  cavity,  relates 
the  average  output  power  to  the  average  power  in  the  electron  beam,  is  given 
by  r)n  =  l/(2nN)  [Ref.  2).  For  the  same  lasing  wavelength,  the  average  power 
in  the  third  harmonic  is  <Pr>3  =  <Pr>1/(3/3).  The  factor  of  Jl  occurs  because 
the  third  harmonic  is  lasing  with  a  value  of  7  that  is  reduced  by  this  factor, 
and  the  average  electron  beam  power  is  linear  in  7.  In  terms  of  the  fiducial 
parameters  (repetition  rate  of  74,375  MHz,  peak  current  of  2  A,  pulse  width  of 
3  ps ,  and  K  =  1.0), 


<Pr>t  =  »;1<Pb>  =  146.76 /Jx,  and 
<pt>3  =  »?3<pb>  =  28.24 /JX, 


where  A  is  in  /im  and  N  =  130.  See  Figure  2. 
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‘  Figure  Captions 

Figure  1.  Typical  cavity  losses  for  multi-layer  dielectric  mirrors  (open 
circles)  and  the  calculated  three-dimensional  power  gain  for  the  fundamental 
(solid  line)  and  third  harmonic  (dashed  line).  Tang's  three-dimensional 
result  for  7  =  350  and  250  are  shown  as  open  diamonds . 

Figure  2.  Calculated  average  output  power  for  the  NBS/FEL.  The  solid  line 
and  the  open  circles  correspond  to  the  maximum  output  power  available  and  the 
actual  output  power  obtained,  using  the  mirrors  and  wavelengths  listed  in 
Table  1.  The  dashed  line  and  the  open  diamonds  are  the  equivalent  result  for 
the  third  harmonic. 
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APPENDIX  K 

Analysis  of  Free  Electron  Laser  Performance  Utilizing  the 
National  Bureau  of  Standards'  CW  Microtron 
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ANALYSIS  OF  FREE  ELECTRON  LASER  PERFORMANCE 

i 

■s 

UTILIZING  THE  NATIONAL  BUREAU  OF  STANDARDS*  CW  MICROTRON 

Cha-Mei  Tang  and  Phillip  Sprangle 

Plasma  Physics  Division,  Naval  Research  Laboratory,  Washington ,  D.C.  20375-5000 

(  Samuel  Penner 

Center  for  Radiation  Research,  National  Bureau  of  Standards f  Gaithersburg,  MD  20899 

Xavier  K.  Maruyama 

Department  of  Physics,  Naval  Postgraduate  School,  Monterey,  CA  939J3 

Abstract 

The  National  Bureau  of  Standards’  (NBS)  CW  racetrack  microtron  (RTM)  will  be 
utilized  as  a  driver  for  a  free  electron  laser  (FEL)  oscillator.  The  NBS  RTM  possesses 
many  exceptional  properties  of  value  for  the  FEL:  i)  CW  operation,  ii)  energy  from  20-185 
MeV,  iii)  small  energy  spread  and  emittance,  iv)  excellent  energy  stability,  and  v)  high 
average  power.  The  1-D  FEL  gain  formula  predicts  that  the  FEL  would  oscillate  at  the 
fundamental  approximately  from  0.25  pm  to  10  pm  when  up-grading  the  peak  current  to 
>  2  A.  In  this  paper,  we  present  3-D  self-consistent  numerical  results  including  several 
realistic  effects,  such  as  emittance,  betatron  oscillations,  diffraction  and  refraction.  The 
results  indicate  that  the  design  value  of  the  transverse  emittance  is  small  enough  that  it 
does  not  degrade  the  FEL  performance  for  intermediate  to  long  wavelengths,  and  only 
slightly  degrades  the  performance  at  the  shortest  wavelength  under  consideration.  Due  to 
the  good  emittance,  the  current  density  is  high  enough  that  focusing,  or  guiding,  begins 
to  manifest  itself  for  wavelengths  >  2.0  pm. 
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Au  FEL  facility  for  applications,  primarily  in  biomedical  and  material  science  research 
as  well  as  for  basic  physics  and  chemistry,  is  to  be  situated  at  the  National  Bureau  of 
Standards1,2.  A  CW  185  MeV  racetrack  microtron  (RTM)S  is  under  construction.  The 
NBS  Accelerator  Laboratory  consists  of  a  series  of  interconnected,  individually  shielded, 

t 

underground  halls.  The  updated  layout  is  indicated  in  Fig.  1.  The  FEL  is  expected  to  be 
operational  by  1990. 

The  major  limitation  of  an  RTM  as  au  FEL  driver  is  that  its  peak  current  capability 
is  lower  than  electron  linacs  which  operate  in  the  same  energy  range.  However,  the  RTM 
is  superior  to  pulsed  linacs  in  energy  spread  and  emittance.  The  RTM  is  comparable  to  a 
storage  ring  in  terms  of  beam  emittance  and  energy  spread,  but  there  is  no  restriction  on 
insertion  length  or  “stay  clear”  aperture.  The  beam  energy  can  be  varied  continuously  over 
a  wide  range  without  significant  loss  of  performance.  In  addition,  microtrons  are  compact 
and  energy  efficient.  Because  of  the  CW  nature  of  the  RTM,  the  generation  of  coherent 
photons  is  not  hindered  by  a  finite  macropulse  length. 

The  original  design  parameters  of  the  NBS  RTM  are  given  in  Refs.  1-3.  The  design 
calculations  indicate  a  longitudinal  emittance  <  30  keV-degrees  and  a  normalized  trans¬ 
verse  emittance1,2  €n  <  10  mm-mrad.  Based  on  recent  measurements  of  the  performance 
of  the  5  MeV  injector  linac,  the  actual  values  of  both  the  longitudinal  and  transverse 
emittance  are  expected  to  be  smaller  than  the  design  values.  The  injector  system  must 
be  upgraded  to  provide  a  peak  current  of  >  2  A  in  3.5  psec  micropulses,  giving  electron 
pulse  length  ~  0.1  cm.  In  order  to  keep  the  average  electron  beam  power  within  the 
capability  of  the  existing  RF  power  system,  the  new  injector  will  fill  only  a  small  fraction  of 
the  RF  buckets  (e.g.,  1/24, 1/120  depending  on  electron  beam  energy).  We  are  proceeding 
with  a  design  of  a  photocathode  injector  system  for  this  upgrade. 

1-D  Free  Electron  Laser  Analysis 

A  first,  order  evaluation  of  the  FEL  performance  of  the  NILS  UTM  can  be  obtained 
from  the  1-D  small  signal  low  gain  formula1.  The  results  indicate  that  sufficient  gain  can 
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be  obtained  at  fundamental  wavelengths  in  the  range  from  10  fim  >  A  >  0.25  y.m.  The 
formula  for  the  electric  field  amplitude  gain  G  in  the  small  signal,  low  gain  regime,  can  be 
written  as 


where  N  is  the  number  of  wiggler  periods,  is  the  initial  relativistic  gamma  factor, 

t 

crR  =  Trr*  is  the  cross-sectional  area  of  the  radiation,  ra  is  the  minimum  1/e  radius  of 
the  Gaussian  radiation  field  amplitude,  I4  =  17  x  10*  A,  I  is  the  current  in  Amperes, 
I\  =  ( \e\Bv,\wl2irm0ci)RMs  is  the  wiggler  parameter,  Bw  is  the  magnetic  field  in  the 
wiggler,  Aw  is  the  wavelength  of  the  wiggler,  Fi  =  J0(b)  —  Ji(b)  for  a  linearly  polarized 
wiggler,  b  =  A'2/2(l+A'2),  v  =  —N A(u>— u\,)/2e  is  the  normalized  frequency  mismatch,  and 
u>„  ~  272c(2ir/A„)/(l+A*2)  is  the  resonant  angular  frequency.  The  function  d/dv{s\n  v/u)2 
has  a  maximum  value  of  0.54  when  u  =  —1.3. 


The  power  gain  can  be  obtained  by 


G„  =  (l  +  G)2-1. 


In  the  low  gain  regime,  Gp  a  2 G.  The  FEL  will  oscillate  when  the  power  gain  is  greater 
than  the  losses  per  pass  in  the  resonator.  The  1-D  gain  formula  is  only  a  rough  estimate. 
It  is  sensitive  to  the  choice  of  filling  factor. 

The  conceptual  design  consists  of  a  linearly  polarized  wiggler  with  a  period  of  A„,  =  2.8 
cm,  and  a  nominal  magnetic  field  amplitude  of  Bwo  =  5400  G.  This  can  be  constructed 
with  a  hybrid  wiggler  design  with  the  gap  separating  the  wiggler  poles  of  g  =  1.0  cm. 
A  wiggler  can  be  constructed  conceptually  in  more  than  one  section,  such  that  a  wiggler 
of  shorter  length  can  also  be  available.  A  shorter  wiggler  and  a  corresponding  vacuum 
chamber  may  be  necessary  for  long  wavelength  operation. 

Figure  2  is  a  plot  of  the  1-D  maximum  small  signal  power  gain  versus  wavelength, 
assuming  a  conservative  peak  current  of  2  A.  The  open  circles  (o)  are  obtained  with 
electron  beam  energies  of  25,  50,  75,  125.  and  175  MeV.  The  solid  curves  are  obtained  for 
the  same  electron  beam  energies,  but  varying  the  wiggler  amplitude  from  0.b£u.„  to  Bwo. 
The  magnetic  field  in  the  wiggler  is  to  be  changed  by  varying  the  gap  between  the  poles 


w 
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from  1.4  cm  to  1.0  cm.  As  the  magnetic  field  decreases,  the  wavelength  of  the  radiation 
decreases,  and  the  gain  is  reduced. 


3-D  Effects  on  the  Gain 

Since  FELs  are  not  actually  1-D,  3-D  effects  will  change  the  gain.  Some  of  the  3- 
D  effects  that  we  will  examine  in  this  paper  are  finite  transverse  emittance,  radiation 
diffraction  and  refraction,  and  some  effects  associated  with  finite-length  electron  pulses. 

We  will  assess  these  three-dimensional  effects  using  a  fully  3-D  self-consistent  computer 
code,  SHERA,  developed  at  the  Naval  Research  Laboratory.  The  formulation  of  the  wave 
equation  is  based  on  the  source  dependent  expansion5  of  the  radiation  field,  and  the  electron 
dynamics*  sure  evaluated  self-consistently.  We  assume  a  waterbag  distribution  in  the  4-D 
transverse  emittance  space,  which  leads  to  a  parabolic  profile  for  the  electron  beam  density. 
Since  the  energy  spread  of  the  NBS  RTM  is  very  small,  it  will  not  be  considered;  and  we 
will  also  not  treat  the  effects  of  pulse  slippage  on  the  gain.  The  radiations  are  taken  to 
have  a  Rayleigh  length  of  175  cm  with  the  minimum  radiation  waist  located  at  the  center 
of  the  wiggler.  Results  for  two  different  operating  regimes  will  be  presented. 

The  effect  of  the  emittance  on  the  performance  of  the  FEL  will  be  more  important 
for  short  wavelength  operations.  Thus,  our  first  example  will  be  for  A  =  0.23  /tm  with 
7o  =  350.  The  pulse  slippage  distance,  NX  =  0.003  cm,  is  much  shorter  than  the  electron 
pulse  length,  The  minimum  1/e  radiation  field  amplitude  waist  is  r„  =  3.57  x  10-2 
cm.  Plots  of  the  power  gain,  Gp,  versus  the  normalized  frequency  mismatch,  i/,  are  shown 
in  Figure  3.  Curve  (a)  gives  the  1-D  estimate  of  the  gain.  Curves  (b),  (c)  and  (d)  are  the 
gains  calculated  from  the  computer  code  for  normalized  transverse  emittance  of  fyv  =5, 
10  and  20  mm  mrad,  respectively.  The  radii  of  the  electron  beams  were  determined  by 
properly  matching  the  beam  into  the  wiggler,  i.e.,  the  radii  of  the  beams  inside  the  wiggler 
is  uniform.  The  matched  beam  radius  condition  is 


where  K#  =  \/2irK/  Xwj  is  the  betatron  wavenumber  for  wiggler  with  parabolic  pole  faces' , 
where  the  focusing  in  both  transverse  directions  is  equal.  The  matched  edge  radius  of  the 


electron  beam  can  be  rewritten  as 


independent  of  the  beam  energy.  The  matched  edge  radii  of  the  electron  beam  are  r„(,  = 
1.77  x  10-2,  2.50  x  1 0-2  and  3.54  x  10-2  cm  for  normalized  edge  transverse  emittances  of 
eN  =  5,  10  and  20  mm  mrad,  respectively.  If  the  emittance  becomes  larger  than  20  mm 
mrad,  the  radius  of  the  electron  beam  will  become  larger  tli&n  the  radiation  spot  size,  and 
the  gain  will  be  substantially  reduced. 

The  effect  of  finite  emittance  on  the  gain  is  negligible  for  A  =  1.25  fim  with  j0  =  150. 
The  pulse  slippage  distance  in  this  case  is  0.016  cm,  and  it  is  still  unimportant.  Figure  4 
shows  plots  of  the  power  gain,  Gpi  versus  the  normalized  frequency  mismatch,  i/,  similar 
to  Fig.  3.  Again,  the  curve  (a)  gives  the  1-D  estimate  of  gain.  Curves  (b),  (c)  and  (d) 
are  the  gains  calculated  from  the  computer  code  for  normalized  emittance  of  =5,  10 
and  20  mm  mrad,  respectively.  Since  the  wavelength  is  longer,  the  minimum  1/c  radiation 
field  amplitude  waist  becomes  r0  =  8.3  X  10“2  cm,  and  the  electron  beam  radii  are  much 
smaller  than  the  radiation  waist.  The  gain  at  A  =  1.25  fim  is  insensitive  to  the  design 
value  of  the  finite  transverse  emittance. 

Figures  3  and  4  also  show  a  shift  of  the  zero  crossing  of  the  gain  curves  obtained  from 
3-D  simulation.  This  shift  comes  from  the  change  in  the  phase  of  the  diffracting  radiation 
field.  It  has  no  real  important  effect  on  the  oscillation  criteria  for  the  examples  under 
consideration. 

Figure  5  plots  the  maximum  3-D  power  gain  versus  wavelength  with  a  peak  current  of 
2  A,  for  normalized  emittances  of  5,  10,  and  20  mm  mrad.  Each  curve  is  obtained  for  the 
identified  electron  beam  energy,  but  varying  the  magnetic  wiggler  amplitude  from  Bwo  to 
0.6 Bwo,  where  the  longer  wavelength  correspond  to  the  larger  magnetic  field.  Normalized 
emittance  is  very  good  in  the  long  wavelength  operating  regime.  In  the  shortest  wavelength 
operating  regime,  the  normalized  emittance  larger  than  5  mm  mrad  should  be  avoided. 

Since  the  current  is  a  function  of  axial  position  in  a  finite  length  electron  pulse,  and 
pulse  slippage  is  unimportant,  the  local  gain  is  a  function  of  the  local  current  in  the  electron 


pulse.  For  tlie  first  example  at  A  =  0.23  pm  with  normalized  transverse  edge  eruittance 
e/v  =  10  mm  mrad,  the  simulations  indicate  that  the  gain  is  proportional  to  the  local 
current,  consistent  with  the  1-D  formula.  For  the  second  example  at  A  =  1.25  pm,  the 
gain  increases  faster  than  the  linear  power  of  the  current.  Figure  6  is  a  plot  of  normalized 
power  gain,  i.e.,  power  gain  from  simulation  divided  by  the  maximum  1-D  power  gaiiT, 
versus  the  frequency  mismatch  at  A  =  1.25  pin  with  normalized  transverse  edge  eruittance 
ejv  =  10  mm  mrad  for  three  different  currents:  (a)  I  =  4.0  A,  (b)  J  =  2.0  A  and  (c)  I  =  0.5 
A.  We  find  that  the  normalized  gain  increases  as  current  increases.  This  can  be  explained 
by  the  self- focusing  or  guiding  phenomenon5,8-11  of  the  FEL.  This  is  most  easily  observed 
in  the  plots  of  the  normalized  1/e  Gaussian  radiation  field  amplitude  radius,  shown  in  Fig. 
7.  For  I  =  0.5  A,  the  radiation  radius  behaves  like  a  free  space  resonator  radiation  field, 
curve  (---).  For  /  =  2  A,  the  radiation  radius  is  less  than  the  free  space  radius  at  the  end 

of  the  wiggler  as  self-focusing  begins  to  show,  curve  ( - ).  If  the  current  can  be  increased 

to  4  A,  the  radiation  becomes  even  more  focused,  curve  ( . ).  The  reason  that  self-focusing 

is  evident  at  such  low  current  is  that  the  emittance  is  very  good  and  current  density  is  high 
throughout  the  interaction  region,  i.e.,  high  beam  brightness  Bn  —  2//(ir2e^)  >  4  x  10® 
A/m2 /rad2,  where  I  >  2  A  and  edge  emittance  cn  =  10  mm  mrad. 

Conclusions 

The  3-D  self-consistent  simulation  results  from  the  computer  code  SHERA  indicate 
that  the  design  value  of  the  transverse  emittance  is  very  good,  so  that  it  does  not  degrade 
the  FEL  performance  for  intermediate  to  long  wavelengths.  For  the  shortest  wavelength 
under  consideration,  emittance  larger  than  5  mm  mrad  should  be  avoided.  Due  to  the 
good  emittance,  the  current  density  is  high  enough  that  focusing,  or  guiding,  begins  to 
manifest  itself  for  wavelengths  >  2.0  pm. 
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Fig.  1.  Updated  configuration  for  accelerator  and  FEL  halls.  The  entire  shielded 
complex  is  located  40  ft  below  ground  level.  Visible  and  infrared  radiation  will  be  directed 
to  a  ground  level  laboratory  (indicated  by  the  dashed  lines)  above  the  UV  laboratory.  _ 

Fig.  2.  Small  signal  poWer  gain  versus  wavelength  based  on  1-D  calculation. 

Fig.  3.  Power  gain,  Gp,  versus  frequency  mismatch,  ts,  at  A  =  0.23  pm  with  -y0  = 
350.  Curve  (a)  is  based  on  1-D  gain  formula.  Curves  (b),  (c)  and  (d)  are  obtained  from 
simulations  with  normalized  transverse  edge  emittances  of  €n  —  5,  10  and  20  mm  mrad, 
respectively. 

Fig.  4.  Power  gain,  Gp,  verus  frequency  mismatch,  t',  at  A  =  1.25  pm  with  = 
150.  Curve  (a)  is  based  on  1-D  gain  formula.  Curves  (b),  (c)  and  (d)  are  obtained  from 
simulations  with  normalized  transverse  edge  emittances  of  =  5,  10  and  20  mm  mrad, 
respectively. 

Fig.  5.  Power  gain  versus  wavelength  based  on  a  fully  3-D  self-consistent  simulation  by 
varying  energy  and  emittance  of  the  electron  beam,  and  the  magnetic  field  of  the  wiggler. 

Fig.  6.  Normalized  gain  versus  frequency  mismatch  for  A  =  1.26  pm  and  edge 
emittance  of  c/v  =  10  mm  mrad.  Curves  (a),  (b)  and  (c)  correspond  to  results  obtained 
with  currents  of  I  =  4.0,  2.0,  and  0.5  A,  respectively. 

Fig.  7.  Normalized  1/e  Gaussian  radiation  field  amplitude  radius  as  a  function  of 
distance  z  in  the  wiggler  with  A  =  1.25  pm  and  edge  emittance  ejv  =  10  mm  mrad  for 
three  different  currents:  (---)/  =  0.5  A,  ( - )  I  =  2  A  and  ( . )  I  =  4  A. 
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Reflection  Matrix  for  Optical  Resonators  In  FEL  Oscillators 


REFLECTION  MATRIX  FOR  OPTICAL  RESONATORS  IN  FEL  OSCILLATORS 

S.  Riyopoulos*),  P.  Sprangle,  C.  M.  Tang  and  A.  Ting^ 

*  Plasaa  Theory  Branch 

Plasma  Physics  Division 
Naval  Research  Laboratory 
Vashington  DC  20375-3000 

The  transformations  of  Gaussian  radiation  beams  caused  by  reflection 

off  mirrors  is  an  important  issue  for  FELs  operating  as  oscillators.  The 

reflected  radiation  from  a  single  incident  Gaussian  mode  vill  contain  other 

modes  due  to  the  finite  mirror  size,  the  deflection  of  the  beam  and 

mismatches  in  the  curvature.  A  method  for  analytic  computation  of  the 

reflection  matrix  is  .developed  by  talcing,  the  convolution  of  the  source 

function  at  the  surface  of  the  mirror  vith  the  paraxial  propagator.  The 

mirror  surface  that  reflects  spherical  incoming  vavefronts  into  spherical 

outgoing  is  found  to  be  a  paraboloid.  Integral  expressions  for  the 

reflection  coefficients  R*”  for  any  incoming  mode  u  into  the  outgoing  u 

pq  mn  pq 

are  obtained  as  functions  of  the  deflection  angle  f,  the  reflected  beam 

spot  si2e  Vq  and  the  mirror  size.  The  coefficient  Rqq  for  the  lovest-to- 

lovest  mode  reflection  is  determined  analytically.  The  spot  size  V  can 

00  ° 

-then  be  selected,  depending  on  the  mirror  size,  to  maximize  Rqq.  The  ratio 
of  the  mirror  size  to  the  spot  size  is  the  dominant  factor  determining  the 
reflection  coefficient.  The  effects  of  deflecting  the  light  beam  enter  as 
small  corrections,  of  first  order  in  the  diffraction  angle  6^  «  1. 

a) 

Science  Applications  Inti.  Corp.,  McLean,  VA  22102 
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I.  INTRODUCTION 


Free  Electron  Lasers  (FEL)  operating  as  oscillators**^  require  the 

Q  IQ 

trapping  of  light  pulses  betveen  systems  of  mirrors  (resonators). 

These  pulses  are  repeatedly  amplified  via  synchronous  interaction  wijh 

electron  pulses  passing  through  the  vigglier.  The  simplest  oscillator 

configuration  is  that  of  an  open  resonator  vith  two  opposed  identical 

mirrors.  The  radiation  vector  potential  for  this  arrangement  is  expressed 

in  terms  of  the  free  space  eigenmodes  A _ (r)  ■  A _ (r)e _  of  the  paraxial 

■n  on  on 

vave  equation11,  where  e^  is  the  polarization  vector  and  . 


A»n<r> 


Wr>W> 

I  zV/2 

II  + 


2R(z) 


(1) 


2  2 

The  exponent  *(r)  a  k[z  +  (x  +y  )/2R(z)]  in  (1)  contains  the  rapidly 

varying  phase  on  the  wavelength  scale  X  -  2n/k.  The  constant  phase 

wavefronts  are  spherical  of  curvature  1/R(z)  >  z/(z  +  b  ).  The  spot  size 

of  the  radiation  envelope  is  V(z)  -  v  (1  *  z2/b2)1/2,  where  w  «  (2b/k)*/2 

is  the  waist,  and  the  distance  z  is  measured  from  the  position  of  the 

waist.  The  amplitude  squared  of  the  mode  drops  by  1/2  over  a  distance 

equal  to  the  Rayleigh  length  b  (also  known  as  confocal  parameter).  Host  of 

the  radiation  is  confined  within  a  cone  parametrized  by  the  diffraction 

1/2 

angle  0^  ■  V/z  *  (X/bx)  .  The  amplitude  profile  umn(rfW)  contains  the 
transverse  spatial  variation,  equivalent  to  a  small  k^t  perpendicular  to 
the  z-direction.  Higher  modes  correspond  to  a  larger  effective  k^, 
responsible  for  the  phase  shift  exp  liS^fz)].  For  a  given  k  the  mode  is 
completely  defined  by  the  two  independent  parameters  R  and  v  (or  any 
combination  of  two  out  of  the  four  quantities  R,  w,  z  and  b). 

The  functions  umn(rjV)  depend  on  the  elected  coordinate  system.  In 


/ 
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rectangular  coordinates  (x,y,z)  they  are  given  by 


Mt^)  •'*f~  •  <2> 

where  Hb>  Hn  are' the  Heraite  polynomials,  »  (2/V2)*/2(it  2m+n  ml  n!)-1^2 

is -the  normalization  factor,  setting  the  total  electromagnetic  flux  carried 

by  the  mode  equal  to  unity,  and  £  (z)  -  (a  +  n  ♦  1)  tan“*(z/b)  is  the  slow 

inn 

phase.  In  cylindrical  coordinates  (r,S,z),  where  tan©  -  x/y,  r  « 
(x2+y2)1/2,  up(r,QjV)  take  the  form 


where  +p  (-p)  signifies  cosine  (sine)  poloidal  dependence,  Lp  are  the 
Laguerre  polynomials,  ap  -  (l/2nV2)1/2  [m!/(m  +  p)!J“1/2  and  Sp(z)  »  (2m  + 

m  in 

p  1)  tan_1(2/b). 

The  electron  beam  is  an  optically  active  medium  that  alters  the 

characteristic  parameters  of  the  radiation  after  each  passage.  During  the 

build-up  period  the  modal  content  and  the  structure  of  the  light  pulses 

inside  the  oscillator  will  change.  A  numerical  method  has  been  developed 

recently  optimizing  the  representation  for  the  amplified  radiation.  In  the 

12  13 

source  dependent  expansion  *  the  waist  size  and  the  curvature  of  the 
elected  modal  basis  is  tailored  according  to  the  driving  source  term.  That 
minimizes  the  number  of  modes  required  to  describe  the  light  beam.  In 
general,  the  curvature  and  waist  size  of  these  modes  does  not  match  the 
curvature  and  waist  of  the  vacuum  eigenmodes  for  the  resonator.  Therefore 
the  transfer  matrix  for  a  given  mirror  must  be  known  for  arbitrary  incoming 
modes.  This  need  stems  from  computational  as  well  as  physical  reasons. 
The  knowledge  of  the  cavity  reflection  matrix  R,  together  with  the  gain 
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matrix  G  through  the  viggler,  is  necessary  in  determining  the  potential  for 
steady  state  operation. 

During  high  pover  operation,  grazing  mirror  incidence  may  be  necessary 
to  avoid  exceeding  the  dielectric  breakdown  limit  for  the  reflection 
surface.  Also,  incase  of  a  high  per-pass  gain  vith  optical  guiding,  the 
waist  for  the  reflected  radiation  would  be  much  larger  than  the  incoming. 

In  two-mirror  resonators  the  reflected  radiation  could  hit  the  wiggler. 
Therefore,'  ring  resonators,  including  three  or  more  mirrors,  must  be 
employed  for  the  deflection  and  recirculation  of  the  light  pulses.  The 
study  of  the  reflection  matrix  must  be  extended  to  include  the  effects  of 
deflecting  the  light  beam,  in  addition  to  finite  mirror  size  and  curvature 
mismatches. 

The  amplitude  profile  of  the  incoming  radiation  will  be  modified  by 

reflection.  A  single  incident  mode  A  fr.)  will,  in  general,  be  partially 

mn  1 

reflected  into  different  modes  Apq(r0)  where  (m,n)  *  (p,q).  This  is  caused 
by  the  deflection  of  the  light  beam,  the  finite  size  of  the  mirror  and  the 
curvature  -  mismatches.  Reflection  into  other  modes  will  affect  the 
interaction  between  the  electron  beam  and  the  radiation  in  a  number  of 
ways.  First  the  rms  radius  of  the  light  beam  vill  change,  affecting  the 
matching  beam  condition.  Second,  the  light  pulse  vill  spread  axially 
because  of  dispersion  among  different  modes,  since  the  phase  velocity 
depends  on  the  modal  number  (m,n).  Third,  different  phase  advances  during 
reflection  among  the  various  modes  may  render  these  modes  out  of  phase  / 
after  a  number  of  bouncings  off  the  resonator.  For  the  above  reasons  the 
fraction  of  radiation  scattered  into  other  modes  will  contribute  to  the 
losses  in  FEL  oscillators. 


V-  ’ 
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II.  OUTLINE  OF  THE  METHOD 

The  geometry  of  the  reflection  is  illustrated  in  Fig.  1.  The 
subscripts  i  and  o  denote  the  coordinate  system  used  for  expressing 
incoming  and  outgoing  modes.  r^  is  defined  with  the  Zj  axis  along  the 
direction  of  incidence  and  r  has  the  z„  axis  along  the  direction  of 
reflection.  The  origins  are  displaced  from  the  mirror  center  by  1^  and  1Q 
respectively,  where  1^  is  the  distance  of  the  minimum  waist  w^  for  the 
incoming  radiation  and  1Q  is  the  distance  of  the  minimum  waist  wq  for  the 
reflected  mode.  A  third  coordinate  system  r$  with  the  origin  at  the  mirror 
center  and  zg  aligned  vith  the  mirror  axis  will  be  useful  in  the 
computations.  Underlined  quantities  r^,  tq  and  rg  stand  for  the  mirror 
surface  coordinates  in  each  reference  frame.  The  transformations  among  the 
.various  frames  are  defined  by- 


x.  «  x  cos  ?  -  z„  sin  *  > 
1  s  l  s  l 


x  ■  x  cos  I  +  zo  sin  |  , 
os  z  s  z 


yi  -  V 


yo  -  V 


Zj  *  zg  cos  2  +  xs  s*n  2  +  ^i* 


z  ■  z  cos  £  -  x„  sin  £  +  1  . 
os  2  s  to 


We  are  interested  in  cases  when  the  reflected  radiation  remains 
focused  along  some  direction  zq  making  an  angle  ♦  with  the  incoming  z^. 
Then  the  reflected  vector  potential  will  also  be  expandable  in  free  space 
eigenmodes  Ap^(ro)  of  the  paraxial  equation  in  the  new  direction.  The 
mirror  surface  generating  focused  reflection  in  the  desired  direction  can 
not  be  arbitrary  but  must  be  appropriately  defined.  The  angle  of 
deflection  $  will  enter  the  equation  defining  the  mirror  surface.  The  other 
surface  parameter,  namely  the  curvature  i/Rm,  is  3  free  parameter.  It 
determines  the  curvature  1/Rq  for  the  outgoing  modes  given  the  curvature 
1/R^  of  the  incoming  modes.  In  case  of  reflection  by  an  arbitrary  mirror 
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surface,  the' scattered  radiation  cannot,  in  general,  be  covered  by  the 
paraxial  modes  that  do  not  fora  a  complete  set  in  three  dimensions. 

Ve  consider  incoming  radiation  of  given  curvature  and  of  arbitrary 
amplitude  profile  A^(r^),  consisting  of  various  modes  (m,n)  with  the  &ame 
R . ( z ) .  If  both  incident  and  reflected  radiation  are  expanded  into 
eigenmodes, 

A‘<ci>  ■  E'n  ‘.nf'i)  • 


A°<ro>  ■  2  Vc«>  • 

p,q 


the  relation  among  the  incident  and  reflected  expansion  coefficients  c1^, 
c°  is  formulated  by 

pq  ' 

c°a  R  c*  ,  (6a) 


o  r.  „mn  i 

cpq  ^  pq  mn  ’ 


where  R  are  the  elements  of  the  reflection  matrix  R. 

pq 

We  examine  the  case  when  the  mirror  dimensions  p  are  much  larger  than 
the  wavelength  X,  X  <<  p  (otherwise  diffraction  rather  than  reflection 
would  prevail).  We  also  assume  that  the  angle  C  subtended  by  the  mirror 
C  «  p/R  ,  where  R  parametrizes  the  radius  of  curvature,  is  small,  of  the  / 
order  of  the  diffraction  angle  0^,  C  -  0^  -  c.  The  j-th  component  of  the 
reflected  vector  potential  at  distance  |rQ  -  rQ|  >>  X  from  the  mirror 
surface  S  is  then  given  by 


A  (v)<V 


eikIV  Ho  I 


IV  V 


A  (v)(Ho)  (n  *  *r)’ 
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In  Eq.  (7)  n • Ar  is  the  obliqueness  factor  vhere  Ar  ■  (r  -r  )/|r  -r  I  and  n 

is  the  normal  unit  vector  to  the  reflecting  surface.  The  surface  element 

ds  is  given  by  ds  «  S{zQ-f (xo,yo) ldxodyQdzo  where  z0  “  f<x0»y0)  is  the 

surface  equation.  Equation  (7)  is  the  convolution  of  a  source  term  As(r  ) 

— o 

at  the  mirror  surfac'e  with  the  propagator  exp(ik|rQ-r0 1)/ |rQ-ro | ,  i.e.,  a 
superposition  of  spherical  waves  originating  at  S.  The  source  term  As(ro) 
is  specified  from  the  incoming  vector  potential  A^(r^)  through  the  boundary 
conditions  and  the  coordinate  transformations  (4).  Ve  will  assume  a 
perfectly  conducting  surface,  vhere  the  incident  and  reflected  fields  are 
related  by 

As  -  -A1  +  2  (n  •  A1)  n  ,  (8a) 


and  n  is  the  normal  unit  vector  to  the  reflecting  surface.  When  the  solid 
angle  subtended  by  the  mirror  is  small,  the  last  term  in  (8a)  is  very  small 
and  the  boundary  condition  becomes, 


.s  i 

(v)  =  (v)  • 


(8b) 


Relation  (8b)  corresponds  to  a  phase  shift  by  it  during  reflection.  It  is 
independent  of  the  wave  polarization,  thus  the  subscript  (v)  is  dropped. 
Cross  polarization  effects,  due  to  the  last  term  in  (8a)  are  discussed  in 
Ref.  10.  Most  of  the  computations  vill  be  performed  on  the  mirror  surface. 
To  simplify  the  notation  from  now  on,  ve  drop  the  bar  (  )  under  the  mirror 
coordinates  r.  Subscripted  quantities  such  as  r . ,  r  ,  r  will  signify  the 
mirror  surface  in  each  reference  frame.  Unsubscripted  coordinates  vill 
denote  the  observation  point  in  the  reflected  radiation  frame  of  reference. 

We  seek  cases  when  the  reflected  radiation  propagates  focused  along 
zQ,  contained  within  a  cross  section  of  dimensions  x,y  <<  z-z  .  The 


r 


a  tj  *s 

expansion  |r-rQ|  »  (z-zq)  (1  ♦  l(x-xQ)*  +  (y-yQ; l/2(z-zo)  j  replaces  the 
full  propagator  inside  (7)  with  the  paraxial  propagator  U  ^  in  that 
direction, 


ro> 


,  -ik(z-z  ) 

ik  e _ 

2n  z-z 

o 


-ik 

e 


(x-xq)2+  (y-y0)2 
2(z-z  ■ 


(9) 


It  is  known  that  the  profile  of  a  given  eigehmode  Amn(xo,yo,zo)  at  zq 
is  generated  by  the  propagator  U^r,^)  acting  on  the  mode  Amn(x,y,0)  at 
z*0.  The  inverse  propagator  U  therefore  reproduces  Amn(x,y,0)  from 
Amn^xo'^o' zo^ '  This  suggests  multiplying  and  dividing  the  integrand  inside 
(7)  by  exp  .[i#(ro)]  /  [1  +  z^/b^J^2,  rec&st^n8  (?)  in  the  form, 


iA<rJ 


rr  / 

A°(r)  =  JJds  e  S(rQ)  e 


i  MO 


o'  o' 


U-k(r'ro)  * 


(10) 


where  the  source  term  S(rQ)  is, 


S(rQ)  =  A1[ri(rQ)]  (n  •  dr) 


(11) 


and  the  phase  d(rQ)  is  given,  in  outgoing  coordinates  rQ,  by 


A(ro)  »  k 


2  2  2  2 

.  ..  ,  xi  <ro>  *  h  (ro>  xo  * 

z.(r  )  +  z  +  -  - 

2Ri(ro)  2Ro(ro) 


(12) 


The  expression  for  d(rQ)  depends  on  the  angle  $  through  the  transformations 
between  the  incoming  and  the  outgoing  coordinates,  Eqs.  (4).  Expression 
(10)  is  the  approximation  of  the  exact  solution  (7.'  to  order 


4-  ‘ 


L  -  7 


2  2  2  2 

l(x-xQ)  +  (y-yQ)  ]/2(z-zo)  -  e  .  It  is  valid  provided  that  the  surface  S 
produces  focused  reflection  along  the  desired  direction.  Otherwise,  the 
paraxial  limit  will  fail  to  encompass  all  the  radiation  contained  in  the 
original  expression  (7). 

The  term  exp[iA(ro>]  is  varying  rapidly,  on  the  scale  of  the 
wavelength  X.  Therefore,  its  convolution  with  the  slowly  varying  source 
term  over  an  arbitrary  surface  will  be  vanishingly  small.  In  general,  this 
corresponds  to  radiation  scattering  where  only  a  small  fraction  of  the 
incoming  radiation  is  reflected  along  the  considered  direction  $.  The 
integral  (10)  will  be  finite  only  vhen  it  is  possible  to  satisfy  the 
condition  A(rQ)  *  constant  over  some  surface  S.  If,  in  addition,  S  is  much 
larger  than  X,  expression  (10)  will  be  finite  only  vithin  a  narrow  angle  5$ 
around  4>.  This  guarantees  that  the  reflected  radiation  remains  focused 
along  that  direction.  Therefore,  a  condition  that  the  exact  reflected 
radiation  (7)  be  fully  covered  by  the  paraxial  limit  (10)  is  that 

A(rQ)  =  const. ,  (13) 


along  the  surface  S.  Accordingly,  the  optical  path  is  the  same  along  the 
various  rays  connecting  an  incoming  wavefront  with  its  mirror  image 
(reflected)  wavefront. 

Requirement  (13)  defines  the  appropriate  mirror  surface  zo=f0(xo>yo?  +) 

for  reflection  in  the  elected  direction.  Expressing  all  quantities  inside 

(13)  in  the  mirror  coordinate  frame,  applying  the  transformations  (4)  and 

2 

using  the  scaling  x  /R  ~  y  /R  -  c  «  1,  z  /R  ~  c  we  obtain  from  (13) 
sm  sm  sm 


z 


s 


1 

2R  cos 
m  L 


2  * 

COS  2  + 


(14a) 


where 
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V- 


(14b) 


1  1  1 

R  "  2R  +  2R, * 
m  o  i 

Equation  (14a)  is  the  analytic  expression  for  a  paraboloid  surface. 

Rm  parametrizes  the  mirror  curvature,  being  positive  or- negative  foT  a 

convex  or  concave  mirror  respectively.  The  surface  is  reflection  symmetric 

with  (zx)g  and  (zy)g  as  the  symmetry  planes;  there  is  no  rotational 

symmetry  around  zg.  Surface  (14a)  can  also  be  approximated,  to  second 
2  2 

order  in  (x  /R  )  ,  (y  /R  )  by  hyperboloids  or  ellipsoids  defined  by 
s  hi  s  in 

(2s  -  E.c°s  2)  -  *s2c°s2  2  -  h*  ■  R„2cos2  2  ’  (15a) 

(zs  *  Vos  I)2  *  X52cos2  2  4  y,2  ■  R.2'052  2  •  <15b> 

All  the  surfaces  become  spherical  in  the  limit  of  perpendicular  incidence  4 

*  0,  and  plane  mirrors  when  R  -»<*>.  Relation  (14b)  defines  the  curvature 

m 

of  the  reflected  modes  from  the  incoming  mode  curvature  and  the  curvature 
of  the  mirror. 

Switching  Eq.  (12)  into  the  mirror-aligned  coordinates  r  through  Eqs. 

s 

(4),  and  using  the  surface  constraints  (14),  it  follows  that 

A(rs)  =  const.  +  o[  kp  ]  , 

m 

where  p  parametrizes  the  mirror  size.  A  more  complicated  surface  equation 

(higher  than  quadratic  in  x,  y,  z)  is  required  to  improve  the  constancy  to 

a  higher  order.  Since  kp  >>  1,  the  approximation  A(r  )  =  const,  is 

satisfactory  for  a  first  order  expansion  of  the  reflection  matrix  in  powers 

of  p/R„,  as  long  as  p/R  <  (kp)-*.  In  case  p/R  >  (kp)-*,  the  slow 

m  m  m 

variation  of  A(x  ,y  )  over  S  must  be  included.  That  introduces  an 
s  s 

additional  contrioution  in  the  reflection  matrix,  known  as  spherical 
aberation. 
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III.  COMPUTATION  OF  THE  REFLECTION  MATRIX 


The  reflected  radiation  is  expressed  by 


A°'r>  ■  I  Id*odyo  o(Vyo>  e 


-‘WV 


U-k(r,ro)  ■ 


where  a(xQ,yo)  *  S [ xq , yQ , ro<xQ , yQ > ] .  Expanding  the  source  e(xo,yo)  in 

terms  of  umn(xo,yo)f 


“<xo’yo)  ■  2  V,<VVV' 

m,  n 


and  exploiting  the  property  of  the  inverse  propagator  U  ^ ,  the  reflected 
vector  potential  A°(r)  at  z  -  0  becomes 


A° ( x , y , 0 )  =  £  Rmn  u  (x,y;W  ),  (18) 

m,n 

where  VQ(z)  =  v  (1  +  z2/bQ2)1/2f  wq  -  (2bQ/k)1/2.  Expression  (18)  is  a 

complete  decomposition  of  the  reflected  radiation  into  paraxial  eigenmodes 

for  incident  radiation  of  arbitrary  profile. 

According  to  the  definition  (6b),  the  R™”  element  of  the  reflection 

matrix  R  is  obtained  from  the  source  term  a  (x  .y  )  inside  (16)  generated 

pq  o  o 

by  a  single  incident  eigenmode  ^pqlr£^r0^'  The  inteSrati°n  *s  performed 
in  the  mirror-aligned  coordinates,  taking  advantage  of  the  existing  / 
symmetries.  The  coordinates  r^  and  ro,  defining  the  incoming  and  outgoing 
wave  functions,  become  explicit  functions  of  x  ,  y  through  the 
transformations  (4).  The  surface  equation  (14a)  is  used  to  express  zg  in 
terms  of  (>:s<ys)-  The  mirror  boundary 


L  -10 


t 


2  2  *  ' 
x_  cos  J  + 
s  z 


(19) 


is  defined  by  the  intersection  of  the  infinite  surface  (14a)  with  the  plane 

2  2 

z  *  const  »  2  p  cos  (^/2)/  R  .  After  the  above  manipulations  the 

S  B 

reflection  matrix  elements  take  the  fora 


lmn 

pq 


dx  „  “.n<W  "p,<V>s> 

s  ys  .  ,  2  ,1/2 

1  ♦  <1 


1  ♦ 


z0 


i  + 


zi  ‘W 


1/2 


i^xs*ys) 


(20) 


x  e 


i6pq<W  -  ‘OvV 


■ 

X 

x  sin  ' 

cos  £  ( 

1  -  jr  sin 

m 

i 

s  2^ 

2  “ 

where 


mn 


<W 


a  u 


mn 


pq 


(x  .yj  » 


“pJWV1  *«]• 


2  2 

Expression  (20)  is  correct  to  order  p  /Rn  . 

Each  representation  of  R  is  tied  to  the  choice  of  the  basis  functions 
umn(r).  In  any  case  R,  as  given  by  (20),  depends  on  four  parameters1^ 

R  a  R($,  a,  y;  l).  (21) 


4>  is  the  reflection  angle  shown  in  Fig.  1.  a  is  the  ratio  of  the  incoming 

to  the  outgoing  spot  size  at  the  mirror,  a  «  V. (l^)/WQ(lo) .  u  =  p/Wq 

parametrizes  the  mirror  size  compared  to  the  radiation  spot  size.  £  « 

V  /R  scales  as  the  diffraction  angle  6,  =  V  /I  multiplied  by  the 

o  m  o  oo 

curvature  mismatch  R  /R  between  the  mirror  and  the  radiation  wavefronts. 

o  m 

The  spot  size  Wq  enters  as  a  free  parameter  because  only  the  curvature  1/RQ 
for  the  reflected  modes  is  specified  by  the  mirror  geometry.  Since  many 


4.' 
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combinations  of  Vq  and  1Q  apply  to  a  given  curvature  according  to  paragraph 

Eq.  (1),  an  additional  selection  rule  for  Vq  is  needed.  Note  that  Vq  does 

not  have  to  match  V^.  This  is  obvious  in  cases  when  the  mirror  size  p  is 

smaller  that  V..  Each  value  of  V  defines  a  complete  set  of  modes  for  the 
1  o 

reflected  radiation, and  an  equivalent  representation  for  R. 

Parameters  ♦,  a,  and  w  can  be  arbitrary.  In  most  cases  of  interest, 
hovever,  l  is  small,  C  «  1,  of  the  same  order  as  the  diffraction  angle  6^. 
The  analytic  computation  of  the  matrix  elements  is  carried  out  by  expanding 
the  integral  (20)  in  powers  of  £, 

R  -  R(0)  +  £  R(l)  ♦  l2  R(2) .  (22) 

The  first  order  expansion  is  performed  in  Ref.  10.  In  this  paper  ve 
review  some  of  the  general  properties  of  R  and  focus  on  the  reflection 
of  the  lowest  mode  Uqq. 


IV.  LIMITING  CASES 

When  the  mirror  radius  tends  to  infinity  (1/Rffl  -♦  0),  or  in  cases 
of  vertical  incidence  on  the  mirror  (♦  »  0),  the  higher  order 
corrections  in  the  reflection  matrix  R  disappear, 

R  -  R(0)  (23) 

in  both  representations.  The  nondiagonal  elements  in  R  stem  from  the 
finite  mirror  size  only.  If,  in  addition,  the  mirror  size  is  very 
large,  u  >>  1,  it  is  appropriate  to  take  VQ  *  as  best  representation 
for  the  reflected  radiation.  The  o  ■  1  limit  yields 


pmn  gmn 

‘pq  *  pq 


(24) 


Thus,  in  case  of  large  curved  mirror  and  vertical  incidence,  or  large 
plane  mirror  and  arbitrary  incidence,  the  reflection  matrix  is  the 
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identity  matrix. 

The  case  a  «  1  is  of  special  interest  £or  arbitrary  angle  of 

deflection  +  and  mirror  curvature  1/R,  as  it  will  be  explained  in  the 

next  section.  For  finite  mirror  size  p  >  VQ,  (  u  >  1),  there  exists  _ 

zeroth  order  nondiagonal  terms  inside  R(0).  Since  R(0)  is  independent 

» 

of  the  angle  of  deflection  +,  the  mirror  size  yields  the  dominant 
contribution  to  the  reflection  into  modes  different  than  the  incoming. 
The  effects  of  the  deflection  of  the  light  beam  enter  to  first  order  in 
K,  R(l),  or  higher. 

As  the  mirror  size  becomes  very  large  and  the  limits  of  integration 

are  extended  to  infinity  the  orthogonality  among  the  various  modes  umn 

becomes  effective.  The  off-diagonal  terms  in  R(0)  become  comparable  to 

2 

the  first  order  terms  R(l)  roughly  when  1/y  -  £  -  6^.  At  the  limit  u  ■+ 

®  all  the  nondiagonal  elements  of  R  are  reduced  to  order  £  or  higher, 

Rpq  -  S  Rpjj(l)  +  oa2),  m  #  p,  n  *  q,  (25a) 

and  the  only  matrix  elements  of  zeroth  order  in  K  are  the  diagonal 

R™  -  R™(0)  ♦  OU2).  (25b) 

mn  mn 

in  both  Hermite  and  Laguerre  represenations.  The  lowest  correction  in 

2 

the  diagonal  elements  is  of  second  order  £  ,  while  the  first  order 
contribution  disappears.  This  is  consistent  with  flux  conservation 
during  reflection  in  case  of  a  large  mirror. 

The  superposition  principle  can  be  used  to  describe  reflection 
from  more  complex  mirror  surfaces.  In  case  of  a  mirror  with  a  hole  the 
surface  integral  (7)  over  S  is  expressed  as  J,  ■  J^l  -  JS2  where  S1  is 
defined  by  the  mirror  surface  including  the  hole  surface,  and  Sj  is  the 
surface  of  the  hole  only.  The  total  reflection  matrix  R  is  given  by  R  = 


*) 

R(Sj)  -  R(S2),  the  difference  in  the  reflection  matrices  associated  with 
mirrors  and  S2  respectively.  The  transmission  matrix  T  through  a 
screen  with  an  aperture  of  area  S  is  given  by  T  <  -  R,  R  being  the 
reflection  matrix  for  a  mirror  matching  the  aperture  S.  The— 

transmission  matrix* for  radiation  diffracted  behind  a  finite  size  mirror 
is  given  by  T'  «  1  -  e*n  R  where  1  is  the  identity  matrix. 


V.  REFLECTION  OF  THE  LOWEST  ORDER  MODE 

The  computation  of  all  the  truncated  integrals  for  finite  mirror 

surface  is  nontrivial.  Most  applications  however  involve  the  (0,0) 

lowest  order  mode  as  the  dominant  mode  in  both  incoming  and  reflected 

radiation/  The  strategy  here  is  to  compute  the  element  R^®  of  the 

reflection  matrix  first.  Then  the  waist  for  the  reflected  modes  V  can 

o 

be  selected  so  that  it  maximizes  R^q.  The  optimum  representation 
condition 


3R 
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00 
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*  0, 


(26) 


puts  the  maximum  amount  of  the  reflected  radiation  in  the  lowest  mode  (a 
different  mode  and  matrix  element  may  be  chosen,  if  desired).  It  is 
pointed  out  that  (26)  does  not  improve  the  properties  of  the  reflected 
radiation.  It  enables  one  to  choose  the  best  representation  in  terms  of 
minimizing  the  coefficients  of  the  undesired  modes  for  the  scattered 
radiation.  Once  Vq  is  fixed  by  (26)  then  the  exact  location  and  size  of 
the  vaist(s)  for  the  reflected  modes  is  determined  by  solving  the  system  of 
equations 
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The  element  is  identical  in  both  representations  since  the  lowest 
order  mode  Uqq  is  the  same  in  rectangular  and  cylindrical  coordinates. 
Performing  the  integration  (20)  yields  R°q  to  first  order  in  t 


e-(l+«2)y2] 


♦  o(c2). 


(28) 


Note  that  the  first  order  term  vanishes  and  the  lowest  correction  is  of 

2 

second  order  in  £  •  The  exact  dependence  on  the  mirror  size  p  is 

parametrized  by  u  ■  p/VQ,  while  a  ■  Vi/Vo  parametrizes  the  ratio  of  the 
incoming  and  scattered  radiation  spot  sizes  at  the  mirror.  The  optimization 
condition  3R°°(0)/3a  -  0  yields,  a2  «  1  +  exp[-(l+a2)p2]  [a4+(2p2+l)a2-l] . 
In  case  the  mirror  cross  section  is  much  larger  than  the  spot  size  of  the 
incoming  mode,  u  >>  1,  a  -»  1  and  the  reflected  spot  size  at  the  mirror 
matches  the  incoming,  Vq  «  W^. 

Large  mirror  size  is  desired  to  maximize  the  total  reflection 

coefficient.  The  reflection  coefficient  is  given  by  «  P0/Pi  uhere  the 

12  i  2 

incoming  flux  is  »  |c  |  -  E  |cx |  and  the  outgoing  flux  is  given  by 


Po  -  |c 


o  1 2 
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In  Fig.  2  ve  plot  for  the  lowest  order  incoming  mode  as  a  function  of 
p'  =  cos(4>/2)  p/ Uq  «  cos($/2)  p.  p'  parametrizes  the  size  of  the  mirror 
surface  projection  into  the  plane  perpendicular  to  the  incoming 
radiation  direction.  The  incoming  radiation  has  a  wavelength  X  =  lp 
(10  cm),  waist  *  10  cm  at  distance  1^  *  2.5x10  cm  from  the  mirror 
and  radius  of  curvature  (at  the  mirror)  R.  «  2.5x10  cm.  The  mirror  has 
a  radius  of  curvature  Rm  =  2.5x10  cm,  yielding  reflected  modes  >  Rq  = 
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2.5x10  (again  1Q  and  vo  depend  on  the  choice  of  Vq).  In  Fig.  3  ve  plot 

00 

the  reflection  coefficients  R  of  the  lovest  order  mode  (0,0)  into  the 

first  25  modes  (p,q)  with  p  <  q  <  5,  as  a  function  of  y' .  The 
deflection  angle  is  90°  and  the  ratio  of  the  spot  sizes  is  l.“ 
Increasing  mirror  size  maximizes  the  diagonal  element  and  minimizes 
scattering  into  other  modes. 

In  Fig.  4  ve  fix  the  mirror  size  y'  *  2  and  the  angle  4>  «  90°  and 
vary  the  spot  size  ratio  a.  The  best  representation,  maximizing 
and  minimizing  R^j  is  obtained  at  a  -  1.  However,  for  small  mirror  y'  « 
0.66,  the  maximum  for  R®®  occurs  at  a  a  0.70  (see  Fig.  5).  Radiation 
reflected  off  mirrors  smaller  than  the  incoming  spot  size  is  best 
described  by  outgoing  modes  of  reduced  spot  size  Vq  <  IT.  Also  note 
from  Fig.  5b  that  for  small  mirror  size  the  total  power  reflected  into 
the  first  25  modes  never  exceeds  BOX  of  the  incoming  flux;  even  with 
many  more  modes  remains  less  than  1.  In  Fig.  6  the  reflection 
coefficients  R®®  are  plotted  as  functions  of  the  angle  of  deflection  4> 
for  fixed  a  »  1,  y'  ■  2.  It  is  seen  that,  for  sufficiently  large 

reflecting  surface  and  good  choice  of  the  spot  size  VQ,  the  reflection 
matrix  is  not  very  sensitive  to  $.  Comparison  of  Figs.  2  and  3  with  the 
rest  of  the  plots  shows  that  the  relative  mirror  size  to  the  radiation 
spot  size  is  the  most  important  parameter  to  determine  the  reflection 
into  other  than  the  incoming  modes. 
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Figure  Captions 


Figure  1  Reflection  geometry. 

Figure  2  Plot  of  the  total  reflection  coefficient  for  the  lowest  order 
mode  as  a  function  of  the  mirror  size  u'  for  ♦  ■  90®  and  VQ  -  V.. 

Figure  3  Plots  of  the  reflection  matrix  elements  .of  the  lowest  order  mode 
(0,0)  into  the  first  25  modes  p  <  q  <  5.  The  magnitude  l^pgi  *s 
plotted  against  the  relative  mirror  size  u'.  The  angle  of 
deflection  ♦  -  90°  and  a  ■  1  (W^  -  Mq). 

Figure  4  Plots  of  (a)  the  reflection  matrix  elements  lRpq  I '  ant*  (*>)  the 

reflection  coefficient  n,,  into  the  first  25  modes,  against  the 

spot  size  ratio  «  for  p'  -  2  and  $  »  90®. 

Figure  5  Same  as  in  Fig.  4  for  p'  -  0.5. 

Figure  6  Plots  of  the  reflection  matrix  elements  lRpq  I  against  the  angle  of 
deflection  4  for  p'  ■  2  and  a  ■  1. 


